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ABSTRACT 


The objective of the STAR project is to test and space 
qualify a new Continuous cycle cryogenic refrigeration 
system for cooling of sensors and electronics which is 
based upon the newly discovered thermoacoustic heat pumping 
effect. The new refrigerator has no sliding seals, a cycle 
frequency of about 300 HZ, and uses acoustic resonance to 
enhance the overall power density and efficiency. This 
thesis 15 2 50559561511234 Specifically With the design and 
testing of the thermodynamic element (or stack), which is 
responsible for the thermo-acoustic power conversion, and 
the testing of binary inert gas mixtures as working fluids. 
Using the refrigerator’s coefficient of performance 
relative to the ideal Carnot coefficient of performance as 
a measure of efficiency, we have achieved a 93% improvement 


over previous designs. 
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1]. =NI RODUCTTON 


ል. BACKGROUND 
tis beans = tne second ain a Series of theses, covering 


different aspects of the Space Thermoacoustic Refrigerator 


AR) project, a great deal of the introductory material 
has already been written. In order for this thesis to be 
memplete, the introductory material that is felt to be 
necessary will be repeated here. With the exception of 
mel —correctiens, Chapter 1I1 and sections À.1 and A.2 of 
Chapter II have been reproduced from the work of 


mreezpatrick (139887. 
ISTOTY 

Thermoacoustics can generally be described as the 
study of the interaction between heat and sound. Scientific 
interest in this interaction is not new. Lord Rayleigh 
(1878 and 1945) discussed various qualitative examples of 
NOUS Eon” O1 Sound by heat. In one of these examples 
he describes investigations ከሃ Sondhauss (1850) of an 
eee er Observed by slasSoloWwers: heating a bulb of glass at 
the end of a hollow rod occasionally led to the generation 
of sound. Faraday (1818) demonstrated that this effect 
occurred With gases other than air. Sondhauss proved that 
the vibration of the glass itself did not generate the 


sound, but he offered no explanation as to what did. In his 


description of Sondhauss’ work, Rayleigh stressed the 
importance of the phase difference between temperature and 
particle motion. Rijke (1859) built an open tube with a 
wire mesh inside. When the mesh was heated, the tube 
produced a sound of impressive intensity. The function 
of this oscillator is slightly different since it regui 
steady gas flow for its operation. A similar SOT 
production effect was observed by Taconis, et al. (1949) in 
hollow tubes immersed in liquid helium (4.2?*K). The Taconis 
oscillations were addressed quantitatively by Yazaki, et 
al. (1980). The work described above deals primarily with 
prime movers. These are devices that convert a temperature 
gradient into sound energy. We call this the classical era 
Of thermoacoustics: 

Interest in thermoacoustics was renewed When sS 
idea of the use of acoustical devices as refrigerators 
appeared. Gifford and Longsworth (1966) described the 
pumping of heat along a surface caused by a periodic change 
in the pressure of the adjoining gas. Such a change Can be 
produced by the oscillations of a sound wave. In their 
experiments Merkli and Thomann (1974) explored the He TT 
and cooling و‎ on the ,Wall of a,.gascfilled 2 mE 
tube. They found that heat uas transported from a region 
near the velocity antinode (or maximum value) CONS 
Sound's standing Wave. to regions near the adjacent pressure 


antinodes.  Thermoacoustic prime movers have also received 
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ሸንያ.ሕፓ O emprton (1976) investigated the excess 
noise of  aeroengines above that predicted by theory. He 
determined that the sound was produced by unsteady heat 
transfer. Each of these - latter three groups of 
experimenters used some theory for comparison, but it was 
mls ۱۱2۵ ۱ ۷2 7 1۷۵۰. None provided the complete quantitative 
theory that would explain their experimental observations. 

The theoretical breakthrough in the thermoacoustics 
field was made by Nicklaus Rott (1969, 1974, 1975 and 
1980). His theory combined the basic principles from 
physics, thermodynamics and acoustics to quantitatively 
describe the effects found in both types of thermoacoustic 
devices: prime movers and heat pumps. Prime movers use a 
temperature gradient to create sound, as in the effects dis- 
cussed by Sondhauss, Rayleigh, Taconis and Kempton. Heat 
mms, Or refrigerators, use the oscillating pressure of a 
sound wave TO produce a temperature difference, as 
described by Gifford and Longsworth and Merkli and Thomann. 
Rott described the effect found by Merkli and Thomann as 
thermoacoustic streaming. It is this effect that makes 
thermoacoustic refrigeration possible. 

Inspired by Ceperley’s (1979) traveling wave Ster- 
۱۳ Cycle heat engine and ic wwe i tative theory, the 
1225م فا لدان‎ 1675277 58612186 Swift, Migliori and Garrett 
(1982, 1983a, 1983b, 1985 and 1986) developed a series of 


thermoacoustic 2 ۱۱ at Los Alamos National 


Laboratory in New Mexico. They investigated the basic 
thermoacoustic effects in both prime movers and 
refrigerators and compared their experimental results to 
Rostdcsesthceorys 
2. Thermodynamics 

We'll digress from history here to explain the 
thermodynamic disti NG ion between prime movers and 
refrigerators. This discussion follows Sears and Salinger 
(1979). AÀ prime mover receives heat at a high temperature, 
1 (the hot reservoir), does work on its surroundings MM 


H 


rejects heat at a lower temperature, Te (the comi 
reservoir), as Shown in Figure l-1l.a. The first law ۵ ther- 
modynamics tells us that the heat flow ) from the hot 
reservoir must equal the work done on the surroundings (W) 
plus the heat flow (Qc) to dXheucodldinessEn E 


Q -W*Q,. (1-1) 


H 
The second law of thermodynamics states that in every 
process the entropy of the universe must either remain 
constant or increase. The entropy change for the iSo ss 
heat transport to/from the reservoirs is defined as the 
heat flow divided by the reservoir temperature. For a prime 
mover the second law can be ritto: 

c or መ“! ( Lae 


The thermal efficiency of the prime mover is defined as the 


Su iS 


ratio of the work output toa 5 ۱ 1 اه‎ 0 ۰ 


"| ىم 0 ) -1 جتن 0 / ( م لاحن © ) دن / انا‎ y1. (LS 
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ከ LAW: Qp = W + Dc‏ ردك 


SECOND LAW: h 2 Oc 
p በ 

EFFICIENCY: 
W ۳ 


7۲ = — ىع‎ |- — 
O. Ta 


FIRST LAW: Qu- W + Qc 


Si 


Ciba ۱ ይስ OF PERFORMANCE : 





O T 
s NEC -:-- 
W 1۳ Te 
5 
Pigure I-1. The Thermodynamics of (a) Prime Movers 


and (b) Refrigerators 


Comparing equations 1-2 and 1-3 gives us a limit on the 
efficiency based on the temperatures of the reservoirs: 
"11-4171 ( (I-4) 


This limit is known as the Carnot efficiency. 


In a refrigerator the process is basically reversed 


(see figure 1፦1.5.). Heat (Qu? is removed from the low 
temperature reservoir, work (W) is received fron ን 
surroundings, and heat Qy? is delivered to the higher 
temperature reservoir. The expression for the first law is 


the same as for the prime mover, Equation I-1. For the 
refrigerator it means the heat flow into the hot reservoir 
must equal the heat flow from the cold reservoir plus the 
work done on the systen. The second law for this case is: 
م255 رتش‎ Or Qu Tyto To: (1739 
The efficiency parameter for the refrigerator is called the 
coefficient of performance (COP). It is defined as the 
ratio of the heat flow from the cold reservoir to the Work 
input to the refrigerator, or: 
of (THT). (TSON 


The limit 15 knoun AS the Carnot coefficient of 


01-0 نا ى‎ 7) $T 


performance. 

The work in this thesis deals solely with refriger- 
ators. Even though prime movers are mentioned from time to 
time, the primary focus of the rest of our discussion uill 


be refrigerators. 


3 Overall Efficiency 


MECO AAN CORO Lera 1988 “3غ 1 عدج مااع ط زا(‎ thermoacoustic 
Work at Los Alamos focused on experimental refrigerating 
engines. The. performance of these engines fell short of 
expectations, leading to simple experiments con basic 
thermoacoustic effects and a proof-of-principle experiment 
on  thermoacoustic refrigeration. Hofler then applied the 
Rott theory to the experimental systems and solved the 
resulting equations numerically. For his doctoral 
dissertation from the University of California, San Diego, 
Hofler designed and constructed a completely functional 
thermoacoustic refrigerator. He also made accurate 
measurements of its thermodynamic efficiency, and used this 
efficiency to make comparisons to the Rott theory. After 
receiving his doctorate, Hofler came to the Naval 
Postgraduate School (NPS) as a post-doctoral fellow, and 
brought his prototype refrigerator with him. 

Them purpose OT this thesis, mY conjunction with the 
work of several other students, is to modify Hofler’s 
PET 1gerator design ima order wo "improve its overall 
efficiency and make it suitable for space cryocooler 
appl ications. The basic design of the thermoacoustic 
refrigerator is shown schematically in Figure I-2. The 
driver (A), “Which produces the sound waves, is coupled to 
the resonator (D) using bolts, via a reducer cone (B) and 


bellows (C). Inside the resonator is a Stack of plastic 
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2 
@ = 3 
1 © 
0 


Figure I-2. Basic Thermoacoustic Refrigerator 
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poke S (E) and their heat exchangers (F). The heat 
exchangers allow heat to be removed from the hot end of the 
ከ aNd applied to the cold end. It is the interaction 
between the sound waves and the plastic plates that 
00 1165 a temperature difference across the plates and 
allows heat to be pumped from the cold heat exchanger to 


the hot heat exchanger. A brief description of efficiency 


would be useful here. Figure I-3 shows a diagram of the 
chain of energy transformations MIG OCCUT IM a 
thermoacoustic refrigerator. There are three conjugate 


variable pairs (V and I, P and U, and T and Qo), and two 
۱۳ ۵ 9 ۲ 1 0 5 telectrical-to-acoustical and acoustical-to- 
۱ (۰ 
We start with electrical input power to the driver, 
لات 1 لنا‎ can be calculated (using r.m.s. values) with: 
عم ی دح‎ (۰ (I-7) 


Where I is the current into the electroacoustic driver, V 


is the voltage across the driver terminals, % is the phase 


Ende “between I and V, and cos?’ is known as the “power 
factor." An electrodynamic driver converts the electrical 
power tO acoustic power (the first transformation). 


Analogous to the previous definition of the electrical 
11: انام‎ poðğer, the acoustic power is given by: 

ዘ=፻ሀር65#, (I-8) 
where P is the acoustic pressure, U is the volumetric 


velocity, Which equals the product of particle velocity and 
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Figure I-3. Chain of Energy Transformations 
in a Thermoacoustic Refrigerator 
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resonator cross-sectional area, and P is the phase angle 
betueen them. The efficiency of this electroacoustic trans- 
formation is given by: 

MEP ec (1-9) 

In the second transformation, the refrigerator con- 

verts the acoustic power to a temperature gradient ( T) and 

Meat flow (Q). As discussed previously, the efficiency for 
this transition is given by the coefficient of performance: 
008-۵ ۰ (1-10) 

The overall thermoacoustic efficiency is therefore the 


product of the efficiencies for the two transformations: 


) (COP) =Q,/P AE 


የጥል” ጅለ elec' 


Improvements in the overall refrigerator efficiency 
can thus be divided into two distinct, though interrelated, 
efforts. The First is the subject of this thesis: 
modifications to improve the coefficient of performance. 
The second effort involves modifications to the driver to 
Optimize the electroacoustic efficiency and is addressed by 
LT Michele Fitzpatrick, U.S. Coast Guard, in her master’s 


thesis. 


B. SPACE THERMOACOUSTIC REFRIGERATOR (STAR) 
1. Motivation 
In addition to improving overall efficiency, our 
experimental refrigerator will be designed and built with 


the intention of flying it on the space shuttle as part of 


|| 


the National Aeronautics and Space Administration’s (NASA) 
Get Away Special (GAS) program. As of 3 February 1988, 
there is a Memorandum of Agreement between the Naval Post- 
graduate School (NPS) and the Air Force which provides 
funding for the space flight and which assigns the NASA 
payload number G-337 to this project. The reason for 
testing the thermoacoustic refrigerator in space is that 
the immediate applications of the STAR are expected to be 
the cooling of electronics, high To superconductors and 
infrared detectors in space. 

There are currently two primary methods for cooling 
satellite systems in space: (1) evaporation of expendable 
Cryogens (liquid helium, nitrogen, ammonia-methane and 
solid hydrogen, etc.), and (2) closed cycle refrigerators 
(Stirling cycle, Vuilleumier cycle, etc.), which involve 
large reciprocating piston masses Operating at lou 
frequency. [Walker (1983) and Smith, et al. (1984)] The 
disadvantages of these current cooling methods are their 
Short lifetimes (expendable cryogens) and high vibration 
levels and low reliability (closed cycle refrigerators). 
The advantages of the thermoacoustic refrigerator 8 
these and other Systems lie in its Simplicity and 
reliability. The STAR has no sliding seals, relatively 
efficient performance, low vibration levels and (hopefully) 


a longer life span. 


The thermoacoustic refrigerator needs to be 
surrounded by a vacuum and insulating material to minimize 
heat loss due to? thermal conduction and radiation. The 
vacuum of space will provide the best insulation to 
determine the refrigerator’s ultimate efficiency. The 
absence of gravity removes the possibility of thermal 
convection as an additional nuisance heat transport effect. 

2. Get Away Special (GAS) Program 

NASA’s GAS program allows small, self-contained pay- 
loads to fly on the Space Shuttle in GAS canisters for 
relatively low cost ($10,000) (Get Away Special Team, 
1984). The GAS can is five cubic feet in volume and can 
house a payload of up to 3200 pounds. Each payload must 
contain its own electrical power, control, data acquisition 
and storage facilities, etc. The Shuttle’s astronauts will 
turn a switch on or off at designated times during the 
flight, but are otherwise not involved with the experiment. 

3. Shared Subsystems from NASA G-313 

Another group of NPS students and faculty (Boyd, et 
31.58 1987) have taken advantage of the GAS program to 
DOI. the» resonant acousticemnodes of the shuttle payload 
bay and the ambient acoustic environment produced as a 
result of main engine and booster operation during launch. 
Their experiment is titled "The Space Shuttle Cargo Bay 
Vibroacoustics Experiment" and is designated by NASA as 


payload G-313. Several subsystems that were developed for 
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NASA 4-313 will be used by "هق شلا‎ 6-337 (STAR). A schematic 
of the STAR in its GAS can is shown in Figure I-4. One of 
the systems borrowed from NASA 6-313 is the controller 
system used to run the experiment and record the data. The 
recorder system consists of the INTEL model BPK 5 
magnetic bubble memory module and an NSC 800 micro- 
processor-based controller. Two other NPS students, LT 
Charles B. Cameron, USN and CPT Ronald Byrnes, USA, will be 
designing the analog electronics and software to integrate 
these systems into the STAR experiment for their master's 
theses. A block diagram of the electronics is shown in 
Figure 1-5. 

Another borrowed system is the power supply, which 
consists of Gates brand lead acid battery cells (five 
ampere-hour, two volts each). These gelled electrolyte 
batteries are ideal for the STAR due to their high power 
density, low cost and the low level of outgassing during 
discharge cycles. NASA G-313 used a one layer battery of 
68 cells providing 680 watt-hours of energy and weighing 
about 80 pounds (including the cells”? support structure). 
We will be using two battery layers with as many cells as 
We can, up to 136. This will give us a minimum of 680 uatt- 
hours and a maximum of 1,360 watt-hours of available 
electrical energy, and a total battery weight between 80 
and 160 pounds. The number of cells ue can use will depend 


on the total weight of the driver-resonator assembly sand 
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Figure I-4. Schematic of the Space Thermoacoustic 
Refrigerator in its Get Away Special Cannister 
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Figure I-5. Block Diagram of the Electronics 
for the Space Thermoacoustic Refrigerator 
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11-5 ጠበ ፣ ከስየየያይመጅጵቹቹቓቻያ equipment (vacuum can, gas reservoir, etc.) 
and electronics, as we have to stay within the 200 pound 
limit. 

The use of the GÁS can imposes certain restric- 
tions. Since we are using batteries to supply the power, 
the refrigerator has to be energy efficient. Also, the GAS 
Can setup requires the STAR to be compact and lightweight. 
Figure e shows a photograph of Hofler’s prototype 
refrigerator. This setup, which includes parts that will be 
unnecessary in space, iS approximately six feet in height. 
In comparison, the maximum payload height for the GAS can 
i5 تبون‎ inches, or less than E. 5 feet. These 
considerations played a major role on the choice of 


equipment for and the design of the STAR. 


C. SCOPE 

This thesis is primarily experimental in nature. The 
purpose being to improve the thermodynamic efficiency of 
Hofler”s prototype refrigerator. Theory is presented only 
as necessary to provide the reader with an understanding of 
the  thermoacoustic refrigeration process and to explain why 
certain modifications were made. For this reason, no direct 
comparison of overall refrigerator performance to theory is 
presented. Hofler (1986) presents a quantitative comparison 
of the refrigerator’s performance to the detailed Rott 


theory. 
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Figure I-6. Photograph of Hofler's ( 1986) 


Prototype Refrigerator 
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Chapter و‎ ۰ ۵9۵ ۳ 3 fairly inedepthediscussion of 
thermoacoustic theory, to provide the above mentioned 
understanding and to illustrate possible areas for the 
improvement of efficiency. This is followed by the details 
of how and why specific modifications, such as gas mixtures 
and stack designs, were made. 

Chapter III. provides all of the information necessary 


for the construction of a thermoacoustic stack. 


Chapter IV. provides a discussion of the overall 
experimental apparatus and procedures. We have basically 
used the . prototype system of Hofler (1986) for our 


experiments. The material in sections A and B of this 
chapter are derived primarily from Hofler’s doctoral 
dissertation (1986), with emphasis on the few areas that we 
have changed. 

Chapter V. presents the methods used for analysis of our 
data and the overall results of our experiments. 

Chapter VI. provides conclusions and recommendations for 


further development. 


Il. THEORY 


A.  THERMOACOUSTIC THEORY 

As stated in Chapter Í., the first tuo sections (A.1. 
and A.2) of this Chapter are essentially reproduced from 
the master’s thesis of LT Michele Fitzpatrick (1988). 

Thermoacoustic theory has been developed in detail by 
Rott (1969, 1974, 1975 and 1980) and adapted to the thermo- 
acoustic refrigerator by Wheatley et al. (1982, ASAS 
[ 9909 1985 and 1986), Wheatley and Cox (1985), Hofler 
(1986) and Sui tana 

The space ی‎ refrigerator (STAR) basically 
consists of an acoustic driver producing sound Waves in a 
resonant tube (see Figure 1-2). This resonator is filled 
With a mixture of helium-xenon gas (412.5 atomic % xenon) 
pressurized to ten atmospheres, and contains a stack of 
plastic plates. This section of the Chapter will present a 
qualitative model for the thermoacoustic heat pumping 
process followed by a quantitative development for heat and 
work flow at a single plate. This is followed by a more 
rigorous development for the heat flow, work flow and 
coefficient of performance for a stack of plates that are 
much Shorter than one quarter wavelength. Section A. 


concludes with a discussion of some of the changes that 
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could improve the coefficient of performance of the 
refrigerator. 
1. Qualitative Picture 

Consider a parcel of gas that moves back and forth 
along one of the plates at the acoustic frequency (see 
mere As it moves, the gas parcel will experience 
changes in temperature and volume. Part of the temperature 
changes come from the adiabatic compression and expansion 
of the gas by the sound pressure, and part as a consequence 
of the local temperature of the plate itself. A temperature 
gradient may develop along the plate as a result of the 
operation of the refrigerator. The temperature and volume 
changes can be described by six separate steps (the changes 
in volume are indicated by changes in the size of the 
square parcel in the figure). 

Assume the plate is at a uniform temperature, T. The 
parcel starts at position 1 (A=Q0) with a temperature of fT. 
The sound wave moves the parcel to the right to position 2. 
The parcel has undergone adiabatic compression and its 
temperature is now T. Since the temperature of the gas is 
higher than that of the plate, heat (Q) will flou from the 
gas parcel to the plate. The parcel’s volume is decreased 
and its temperature lowered to q”. This heat flow also 
causes the plate’s temperature to increase at the position 
x. Parcel position 3 is actually in the same location as 


position 2, but is displaced vertically in the figure for 


ሯቹ) 





COLD END HOT END 


EA TE 


Gas Parcel Diagram (Note: ۲۱5 8 
in the Vertical Direction 15 Shown for 
011807777095595 ግ ን Horizontal 
Direction e. 
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Figure ] 


0001000 Ihe parcel now moves to position 4 and expands 
adiabatically to a new temperature of T . A repeat of this 
Process puts the parcel in position 5 with a temperature of 


1 . Nou the temperature of the gas is louer than that of 


the plate and heat (Q’) flows from the plate to the gas 


parcel, "expanding the parcel and raising its temperature to 
T- . This heat flow causes the plate’s temperature to 
decrease at position qe. The parcel moves to the right, 


again under adiabatic compression, and we are back where we 
started, in position 1 with temperature T (adapted from 
Wheatley, et al., 1985). 

This gas parcel cycle is repeated at the resonator’s 
operating frequency. For the STAR this frequency is between 
approximately 250 and 600 cycles per second, depending on 
the gas types and mixture (specifically the speed of sound 
in the gas) and the length of the resonator. 

Notice that the temperature of the gas parcel at 
position X=0 is different depending on which direction the 
gas parcel is moving. It is this phase shift in temperature 
relative to motion that produces the thermoacoustic effect, 
as we Will show in the calculations that follow. 

It is also important to note that in order for heat 


to flow between the gas parcel and the plate, the parcel 


must be vertically located within about one thermal 
penetration depth of the plate. The thermal penetration 
depth (64) is the distance that heat can diffuse through 


e 


the fluid during a time lw, where w is the acoustic 
angular frequency. 

Wheatley, m etann (1986) describe thermoacoustic 
engines as consisting of long trains of these gas parcels, 
all about a thermal penetration depth from the plate. The 
parcels drau heat from the plate at one extreme of their 
oscillatory motion and deposit heat at the other extreme. 
Adjacent heat flows cancel, except at the ends of the 
plate. The net result is that an amount of heat MINS 
passed from one end of the plate to the other. It should be 
noted that the above discussion assumes zero viscosity. 

2. Simgde Plate 

owift (1989) illustrated the basic principles of the 
thermoacoustic engine by examining a simple example of a 
single plate in a gas-filled tube such as we described 
above. In this example the acoustic and thermodynamic 
effects are nearly distinct. In the absence of the plastic 
plate, the plane standing wave sustained in the tube has 
well known acoustic properties. It’s pressure and x- 
direction velocity (along the length of the tube) are given 


as (using peak values): 


p,-P,Sin(x/x) (31:8 
and << )1 7۵۵۵ ٩ (115% 
where P is the pressure amplitude at the pressure 


A 


antinode, *=A/2 T=a/w is the reduced wavelength, A is the 


wavelength, a is the speed of sound in the gas, w is the 


24 


angular frequency and pa is the mean density of the gas. 
The subscript 1 indicates the first order time oscillation, 
and the subscript m indicates the mean value. The total 
pressure can be given as: 


p^ (T3) 


p=p, +p,e” 
A similar expression may be used to represent the total 
ክው ከ5 ስ መ ህ1በመ “Swift, we will assume that the first 
order in the acoustic amplitude is sufficient for all 
acoustic and thermodynamic variables (pressure, velocity, 
temperature, density and entropy). 
The: sound wave is adiabatic and has an oscillatory 
temperature given by: 

e (II-4) 
where A e wobarTc thermal. expansion 
coefficient and So is the isobaric (constant pressure) heat 
capacity per unit mass. Note that T, and p, are in phase. 
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For ideal gases, thermodynamics tells us that: 


=(፻፦ ፦ 
an aS ( DT TP, ۱ 
Where NM is the ratio of isobaric t0 isochoric 
specific heats. Y equals 5/3 for monatomic gases and is 
smaller but still greater than one for other gases. 


Combining Equations II-4 and II-5 gives us: 
፡1/፤እ2(፡ሻ-1)/#3ጅ1/“ጅእ- Gils) 
The introduction of the plastic plate into the 


Standing wave modifies the original, unperturbed tempera- 


ture oscillations. This modification is due to the heat 


ሯን 


flow betueen the gas and the plate, as described 
previously. The temperature is modified in both magnitude 


and phase for gas about a thermal penetration depth away 


from the plate. According to Swift, this results in two 
important effects: 1) a time averaged heat flow near the 
surface of the plate, along the direction of acoustic 
vibration, and 2) the generation or absorption of real 
acoustic power near the surface of the plate. In 


refrigerators the effect is an absorption of acoustic 
power. In prime movers the acoustic power is generated by 
the temperature gradient present in the plate. 

There are several assumptions we will make to 
simplify calculations, again following Swift. We assume 
that the plate is short enough compared to a reduced 
Wavelength (ልጁ«ኋን and far enough from both velocity and 
pressure nodes that P4 and u, can be considered uniform 
over the entire plate. We assume the gas has zero 
ViSCOSity, 50 that u, does not depend on y (ue had already 


1 


made this assumption ከሃ looking at 9 1 in the x-direc am 
only). We assume that the plate has a large enough heat 
capacity per unit area that its temperature does not change 
appreciably at the acoustic frequency. We assume that the 
plate has a mean temperature gradient in the x-direction 


VT ۰ Finally we neglect the plates and gas’ thermal 


conductivity in the x-direction. 
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Applying these assumptions, we see that the mean 
fluid temperature CT af)? is the same as that of the plate. 
Swift calculated the oscillating fluid temperature using 
the general equation of heat transfer. He kept only first- 
order terms, neglected thermal conduction along x and 
applied the boundary condition, T,(0)-0, imposed by the 
plate. The resulting equation is: 

A A መሮ ASH 76.24 . (II-7) 
To interpret this equation, we look at it in the limit that 


the gas is far enough from the plate .(y»$,) to make 


negligible thermal contact with the plate. This gives: 


ነሽ. ክባ SIA, ued. (11-8) 
The first term in Equation II-8 is the temperature 
oscillation due to the adiabatic expansions and 
compressions of the gas (see Equation II-4). The second 


term comes from the mean temperature gradient in the gas. 
As the gas oscillates along x with displacement amplitude 
/w, the temperature at a given point in space oscillates 
by an amount VIR die e even if the temperature of a given 
piece of fluid remains constant. The actual temperature 
oscillations are a linear superposition of these tuo 
effects. 

The y dependent part of Equation II-7 is complex. It 
approaches 1 for y»?é,, as given previously. It approaches 
zero for y*é,, Where the plate imposes the condition 1,=0. 
Most importantly, for yê, its magnitude is still on the 


En 


order of 1, but it has a substantial imaginary part. This 
phase shift in the oscillating temperature of the standing 
Wave at y*é,, due to the thermal presence of the plate, is 
an important result because it leads directly to the time 
averaged heat flow in the x-direction. This is the same 
phase shift that was emphasized in the description of the 
gas parcel motion. 
a. Heat Flow 
Swift argues that since we are neglecting ordi- 
nary thermal conductivity in the x direction, the only way 
heat can be transported along x is by the hydrodynamic 
transport of entropy, carried by the oscillatory velocity 
uj: 
q) 5 n g31947 6172) f, Sp ImCT 11፡1 (PSA 


The line above the quantity s indicates the time-average 


m 
of the first order entropy and velocity product. The heat 
flow is a second order quantity, signified by the subscript 
2. 

The total heat flou 92 along the plate, in the x 
direction, is found by integrating 92 over the y-Z plane: 

"m. 
zm 9247 ) [ 1-10 ( 

where TÍ is the perimeter of the plate in the y-z plane. 


Substituting for T, and performing the integration gives: 


1 
0 -<و‎ )1/4(118 (T زلار ۲( و‎ )۲-1( (11-11) 
where پر‎ is the thermodynamically active area in a plane 


perpendicular to the longitudinal acoustic motion, Tm? is 
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the heat parameter of the gas, and Fl is defined as the 
ratio of actual temperature gradient to the critical 
temperature gradient OTT). The critical temperature 
gradient occurs uhen the temperature change along the plate 
just matches the temperature change due to the adiabatic 
compression of the gas, and no heat flows between the gas 
and the plate. This is the boundary between the 
refrigerator and prime mover functions of thermoacoustic 
engines. 

Equation II-11 shows that when [f<l, heat flows 
up the temperature gradient from cold to hot and work 
(acoustic power) is absorbed, as for a refrigerator. When 
lr] there is no heat flow. When r>1, heat flows down the 
temperature gradient from hot to cold and acoustic power is 
produced, as for a prime mover. 

Note that the total heat flow is proportional to 
the area ff&, and to T. ^ (which equals 1 for ideal gases). 


It is also proportional to the product pu and so equals 


1" 
zero if the plate is at either a pressure or velocity node 
of the standing wave. The maximum value of the product 
occurs halfway between the nodes. Finally, the heat flow is 
proportional to the temperature gradient factor rnr-1. For 


vales T r-1 is greater than zero and the heat flow is 


Crit’ 


r-1 is less 


toward the pressure node, while for VT_<VT 5 
m Crit 


than zero and the heat flow is away from the pressure node. 


If suitable heat exchangers, at Tu and Te: are installed at 
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the ends of the plate (with Ty To= T OX) this heat flow 
carries heat from one exchanger to the other. 

The heat flow is small under ordinary 
circumstances. However, in closed resonators it is possible 
to achieve sound amplitudes many orders of magnitude higher 
than those of ordinary conversation. Since 92 3= 
proportional to Pats and since in practical acoustic 
engines the entire cross section of the standing wave is 
filled with plates (spaced roughly 4$, apart), very high 
heat flows and/or a large Ty -Te may be achieved. 

b. Work Flou 

The work flow (i.e. acoustic power) is given by 
the work per cycle times the rate at which that work occurs 
(the acoustic frequency f). From thermodynamics, the 
average acoustic power produced per unit volume is: 


u =- (w/g YIP, ¢y=-( 1/2) wep ImCT ۰ (11-12) 


1 

The gas about a thermal penetration depth auay 
from the plate "breathes", because of thermal expansion and 
contraction, with the right time phase with respect to 
oscillating pressure to do (or absorb) net work. This is 
exactly the same gas that we have seen is responsible for 
the heat flow. Gas elsewhere is ineffective in doing (or 
absorbing) work. The density oscillations for y«s, and for 


y>6; are in phase with the pressure oscillations, and hence 


do (or absorb) no net work. 
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The total acoustic power produced, Wo. is found 


by integrating ሠጋ Over all space, as with the heat flou: 


e 


U ) 1/4( 1۲8, )۷-1( ) ۳۱۶/۵۰۵۱۲۰۰, (11-13) 


2 
Where all the terms have been defined previously. 

The acoustic power is proportional to the volume 
ሸለመ of a fluid that is about a thermal penetration depth 
from the plate. It is proportional to Pis and so is 
quadratic in the acoustic amplitude (as was the heat flow) 
and vanishes at pressure nodes. Finally ዘጋ is proportional 
to f-1, the same temperature gradient factor as appeared in 


0 (Equation II-11). When INVE r-1-0, and there are 


CENG 
no temperature oscillations in the fluid, other than those 
due to adiabatic processes, and no acoustic power is 
absorbed or generated. For UD Hes r-1»0 and acoustic 
power is produced near the plate. Whether this power 
increases the amplitude of the standing wave, is radiated 
away to infinity, is simply absorbed, or flows through an 
acoustic-to-electric transducer to generate electric power 
depends on details of the resonator, not on the plate 
itself or the standing wave near the plate. For VAT دي‎ 
r-1«0 and acoustic power is absorbed near the plate. For a 
tube without plates, which has a diameter less than the 
wavelength, at constant temperature (VT-O), this work flou 


is responsible for the ordinary thermal attenuation of a 


sound wave (Kinsler, et al., 1982). 


E 


C. Efficiency 
We combine Equations II-11 and II-13 to get the 
efficiency of the plate with no viscous losses or 
longitudinal conduction losses: 
= : i = = - z 
n 4/0, pe: 1)/(፲ 8) 3(ሠለኗጅ))/(የ;ኣ uy? ( 1 1-14) 


Since u,-ugsin(x/2»X), and p,7 Qau cos(x/x), we get: 


m 0 
n=[(7-1)/(Tmp)J3(4x/(1taníx/1))). (11-15) 

For x<7, this reduces to: 
۷ <] ) ۲-1۸ ) ۲11۵ ۶ ( [ ) ۵2۸/26۱ ب‎ ) ] [1 - 16 


and mr, (11-17) 


="carnot 
This efficiency approaches the Carnot efficiency 
as the power output and heat transfer rates approach zero. 
We can make a Similar calculation for the 
refrigerator mode of operation, Where the relevant 
efficiency is the coefficient of performance (COP). Given 


that 008-0۷ we find that: 


2 


COPS dms (113189 


3. Short Engine 


Having developed the simplified case of a single 
plate we now turn our attention to a more realistic model, 
With an entire stack of plates. Since Swift (1989) has 
incorporated the thermoacoustic theory of Rott (1969, 1973, 
1974, i975 and 1980), we will once again follow his 
development closely. The assumptions that the length of the 
Stack is much shorter than a wavelength and that the 


temperature spanned is much less than the absolute 
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miperarcures mill be retained. However, this more realistic 
model includes viscosity, longitudinal thermal conduction, 
finite (instead of infinite) plate heat capacity and a 
ED 01 Pilates (instead “Or just one). The goal in our 
HEvelopment is to include these effects, which produce many 
Complications, in such a manner as to maintain a connection 
to the IS po car more “easily understood, 
development of the previous section. 

Swift defines a co-ordinate system where the x-axis 
00 50112 the direction of acoustic vibration (longitudinal 
اا‎ the SEED and the y-axis iS perpendicular to the planes 
wm The plates (using a parallel plate geometry) With the y- 
origin centered between two plates. This gives a plate 
spacing of 2Y y and a plate thickness of 21. He begins with 


۳۳۱۳ 2۱۱2 ۲ 1 ۳ Of motion and a continuity equation for the 


۲۲ 0 : 
A A “ከተህ ጭው ۴ پ ۷ 0۳۰ رت‎ ARTS y 
and ۲۳۵ ۳ زا‎ 2)97(0 LE 400 
Ante “e. is density, V is velocity, p is pressure, v is 


kinematic viscosity (Swift uses P for this parameter) and 3 
دم‎ ۱ VISCOSITY. He imposes the boundary condition that 
ህ=ኃ at the fluid-plate interface. He then uses the equation 


for heat flow ina fluid and solid: 


PTCIs/dt+V-7sI]=V°(AVT)+( quadratic velocity terms), ESI 
and ኢር (251 ASADA ET, ANE) 
5 5 5 5 5 
Where A iS thermal conductivity, =5 iS entropy per unit 


E 


mass, T is absolute temperature and the subscript s refers 
to the solid. The applicable boundary conditions are that 
baile and AdT /d y at the fluid-plate interface. 
Equations for each of the needed parameters are then 


linearized, assuming that they all oscillate at the angular 


frequency w: 


۲-۳ + 1. (11-23) 
+ )و‎ t, )11-24( 
ا ل مل وين دن‎ (II-25) 
T-T (+۲ ^5, ) 
Tain ^ (11 
and s-s, GO *s, RE (LIS 
where the subscript um denotes a mean value and the 


Subscript 1 refers to a first-order quantity. Swift then 
nanipulates these equations, by integrating With respect to 
Vs into the form of a wave equation for P, x) and an 
equation for energy flux H, (which includes both heat and 


2 


work since H=Q+W) along x. The wave equation is: (11-29) 


| E 
a O CE do, dp, a: Lo dT dp, 
+ SSS) سس‎ ۳۳ ee ee 
Lt 5. ል የ dX| wi (1=0)(61+£_) dx و‎ 
- 5 m - 5 


and the equation for energy flux is: 





2 
x የ EL 
"76 |21 ماس‎ f, | 


| سس 
2w? | dx ۱ ) 1+ : ۱ )1+ (۱‏ 
m |! S |‏ 
ل TTy,C dT dp, dp ርይ ያ ስበር ር‏ 
በ ር a Sein O‏ 
CE GL r)‏ مدال Zw ۳ dx- d%‏ 
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dT. 
MEO 


We have used ل‎ as Prandtl number, 


fi=tanhl (1+1)y¿/53/L (1+1)y¿/% J, 11) 

ሂኒ, =ቲ8ከከቪ (1+1)ሃ65/*, 3]/1(1+1)ሃ6/», 3) 9 

2 | s ፡ | 4 -33 

and ነ በ 4 tanhl (1+i)y,1 ፡ ¢ 2 ۸ ከ ከ 1፡22) 


where å iS a Specific penetration depth denoted by the 
BMSpBODFriate subscript, 4 for thermal penetration depth in 
the fluid, s for thermal penetration depth in the solid and 
DU Or "viscous penetration depth in the fluid (which is 
roughly the. distance that shear momentum diffuses in time 
۱۳ Swift then applies some simplifying assumptions to 
muse eguations that can be related to the results of 
Perron 2. A boundary layer" approximation (i.e., Yg K 
and o) sets the hyperbolic tangents to unity. Since in 
McCain eines the "plate Spacing is 7 to 4 thermal 
penetration depths (or greater) and the hyperbolic tangent 


function is within a few percent of unity at Yg72 *« و‎ nis 


15 a reasonable approximation. A SHOT SURER. 
approximation (i.e., Bx%A) yields relations for Py and uy 
that SL ICO. those ~ot Section 2, With the 
፣ ۱ 1 1 Of "a ‹1+1/ሃሪ፣ termein Uy ፡)ከህ፡። ን የቢ ስበ for the 


ae = Verne nignher i1lsiae the Stack. “This approximation 


also gives us the fact that م‎ u, ana all thermophysical 


1 
Potes sare independent Of x Within the stack. Finally, 


IE nu 71 5050551 0 assumption iS made (although viscosity 


ES 


will be added back in later) thereby setting v and 1j ۲ 
to zero. Using these approximations Equations I1-29, I1-30 


and II-33 can be reduced to: 








o at d[ 1 dp EREET 
و‎ + -|፦ سس ات‎ )۲-1( 2. (11-34) 
w* dx m dx (Itid cli ا‎ 
፲ ጾ cun 
m Swe 8 m | 
۲] پر - -2< پر‎ Py "1 a A E (II-35) 
ES 1+ و‎ መ 
5 
| nC p? 
and 5 ; ( I 173659 
S E 
¡RA 


where all terms have previously been defined. It can be 
Seen that the first term in Equation II-35 is identical to 
the hydrodynamic heat flux 92 derived for the single plate, 
except for the 1/)1+ term. This term accounts for the 
finite heat capacity of the plate and modifies the ۱ 
condition for T4 at the fluid-plate interface to a non-zero 


value, thereby reducing the heat flux. For the Short TS 


Wo is small ‹#ኃ/ሟጋ= d ولي ال‎ SO heat flux and 


energy flux are approximately equal. The second ۶ በክ 


carnot 


Equation II-35 accounts for ordinary conduction oi wie 
down the temperature gradient by the fluid and 3 ۴ 

The next step is to derive an equation ممم‎ 
acoustic power. Since net acoustic power cannot enter or 
leave the engine in the y-direction the acoustic 
generated or absorbed must involve the difference in 
average acoustic intensity pu, between the two ends of the 


stack. Using this and the fact har OS /dx==iwQ Uy Swift 


1 
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۱-۱۱۱ >23 11 11-34 to arrive at a relation for 
acoustic power: 
1 (2—1)w 


-Më Èx 
4 $ 2 ۶ )14+ ۶ ) 
m D 


W 





>= n ከከ |] 


ከከ  < the same as Equation li-13, with the exception of 
3 ۰ 1 ۳ which Was discussed above. From these 
results ፤ 1 is clear that, under these simpliftying 
0000 10100135 the short engine performs in a manner similar 
5 ۲ predicted by the single plate results. 

We now follow Swift in repeating this derivation, 
put include nO ESI 7150:0551 1 inclusion of viscosity 
Beads. Some complexity to the situation. Under these circum- 
stances u becomes a function of y and the concept of 


1 


۱3 )ا‎ 031 temperature gradient becomes les 


n 


; Well defined. 


L 


Since Ti IE Enciso a resultar fashion, there 
۲ በወ natural definition of الح‎ Houever, ue uill use 
S 3 
AS j 3u i uev 7 50 
Edid Swift) | ከው iu and ci d PIER Neu D 


Where the superscript s defines the value for the standing 
Ehe wewadditaon —of^viscosity 100 11 1 ۵5 the relation for 
dp,/dx CO: 

an” 
A _ fn 1 


= ج س‎ መ) WMA rS 
dx Iu ld.) 


dp, 


Euch ispcomplex even though <u,”> is real. 
Mie تلد‎ ito —sauations 11-29 and 11-30 and 
applying Be boundary layer approximation yields the 


following relations: 


5 


S S K ۱ 
z +ሆ “+ሆተሆሂ 
1 ሚሪ <u, > 1 d a 

















| و‎ 2 - 5 ۳ — 
+ Cl+e 3 (140) (1-5, /yy +5 /Y gy) * 1+2 
ds 
= Cy gat la )--- ) [1 1 2۳ 
5 
dx 
and 
p ai  dfi-f, dp ۲ - 6 P 
Aa فك‎ “ሠ | A A a. 
we dx ls. Beiro oM MD DNE UE 


For the zero viscosity case, f=@, these relations reduce to 
Equations 11-34 2ه‎ 11-35 
Work flow is calculated, using the short omi 


approximation, from the uave equation. This gives: 


ም )۷ س(1-‎ ) )፡ r 
4 e ( 1+7 ( ) 1- وت خيرم لا بربيزوة.‎ EO 
1 م ند‎ ረኳ SL 
= a, 6% A (11-41) 
4 1-5, /ሃሬ? 5, ፤/2ሃሬ፤ 


Which reduces to Equation II-37 uhen 5-20. 

From these relations it is clear that viscosa 
makes the situation much more complex. Since the goal MN 
expression for the coefficient of performance, items 
necessary to determine 32 /W.. We know that H=Q,+W,, such 


2 


that CORS 3.. 4 لنت‎ #2)/8,,=(82/8#3)- dus In order to maintain 


the essential effects of viscosity, while reducing the 


complexity of the resulting equations, we will include 


i 


“4 
us 


viscosity only to the lowest order present, that iS 5 
Additionally we will assume that Yo” $k» which is the case 


for most practical engines as previously discussed ss 
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7 
eves US that oa ገ መ Applying these approximations 


ህአ መ CuHations 1/|-40 and II-41 yields: 


. 1. ፐ #1 ‹81.> 4T, 
ሸኃ=- -ሾ5« ሙ-፦፦-፦መ=ጭ( የ>1)-ሸ(ሃልእ+1 እ.) መ> ; a 
2 e ۱ ^ dx 
and 
pU: )۷ س(1-‎ ) )፡ 1 ۲ 
SA == (MSG xp, KU, pu (II-43) 
4 e 22) 1+ 2 ( 4 
m ጋ 


Now we can take the ratio of these two relations and 
get a useable, understandable relation for the coefficient 
of performance, Which iS What we are trying to improve 


experimentally. The result is: 


Sew insyteep <u, >> ¢r=2) dT. | 
m- il 1 . ሽ 
i. AS ک(‎ | 
Gu o CI) dx | 
COP- = A 
(1/48 ,8x (Y D)w(p )፡(የ-1) 1 ` | 
aa a a a >? 
Pa? Cite.) 4 
The above relations are approximate expressions, which 
should apply to realistic situations, With the only 
remaining assumptions being 1) the SOT stack 
approximation, 2) the boundary layer approximation, and 3) 


a Small Prandtl] number approximation. We have previously 
explained that the boundary layer approximation is 
reasonable for the plate Separations of our stacks. The 
short Stack approximation is ESO mot particularly 
!ا‎ 1 1 172 101 1 55115153167 engines. With a Prandtl number of 
approximately 0.027 for helium-xenon ue have a relatively 


small Prandtl number, so this approximation is also very 


E 


reasonable. In his comparison to theory, Hofler ۲ 
shows that in relatively "long" acoustic engines Equations 
11-29 and 11-30 predict the performance within 29 (tame 
percent. 
4. Design Considerations 

The purpose of this section is to discuss the 
primary modification that we have made based on theory, 
that is, the reduction of the Prandtl number of the system. 
To evaluate how this affects performance we will first look 
at Equation 11-44. We know, from experience, that the 
second term in the numerator is small compared to the first 
term. We also know that * is, in general, small comm 


S 


to unity, so we will set ¿_=0 for evaluation purposes. The 


Product Tae is known to equal unity for ideal gases. 
Applying these conditions, With Swift’s Bs o^ and Yo 4 
we can reduce Equation II-44 to: )[1 1-4 


ኣኣ 
CLA My ¿Py CO 2(የ-1)3/(1-ሀ )3-(ሻ(ሃልእ+13 AT /dx] 
,عمج‎ eee 


CLHMy Ax (¥-1)p «።”“የአ‹ክ) ፦! 
It ais clear from this thal as p Seen 1- م‎ is 
increased resulting in an overall decrease in the numerator 
of our relation, and ው is reduced resulting in a decrease 
in the denominator of our relation. These two effects when 
evaluated fully, show that COP will decrease with reduced 
Prandtl number. However, we know based on our experiments 


that this cannot be true for our system. This effect that 


is apparently the opposite of what we expect, and know to 
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due to some (or one) of the‏ ع:01 565181 15 mue ١‏ ع5 
approximations that we have made being invalid for our‏ 
system. AS ve are about to show, the relations given in‏ 
Equations IM 9 and ( ۲ 5 ۳ yield the proper result,‏ 
inclusion of viscosity only to the‏ کت NI‏ 
lowest order does not provide adequately for the effects of‏ 
viscosity in our system.‏ 


Looking at Equations 11-39 and [I-41 we see that 


they are fairly complex. We will use the simplifications 
oie =07 TEM and ۰۵۰۶/۳: and we will consider the 
3 n n n 


thermal conduction term in the heat flow to be negligible. 


We also know that é,=8,/c0°, and we will use 51 =P con kkmaNd 


0 


> 2۵) 17۷۱ ) 8/۲ sink tithe: Simplify our relations. 


This allows us to re-write Equations Il-39 and II-41 as: 


y, 
acoskxsinkx 1 + 5 + 7 


AH fp 
 .. — —— | ARAS uL 2) 
Po لات ی ا‎ 0150 1+7“ 
(II-46) and 
4ዝ Yp ۰ 
- 2=(Y-1)wéxcos kx ا تت اس‎ 1 
Pg Ne (1+« )(1>5,/ሃኔ፥5( 5» “ሃዕ) ' 


ፍ 
4 "a a p 
vaT Sin kx 


- مس‎ ው (11-47) 
L1-$6y7ygt S CS, yo) * 1 
Now ue can combine Equations II-46 and II-47 to get a 


| 015 101 for COP as follows: 


COSkxsinkx 
Cole 





Excel. ( 11-48) 
5፳(14ሀሆ) (1-5, /ሦሃ5፥።=( 5,»/ሃሬዕ) ( 


4| 


where the brackets t---3 represent the following relation: 
4 2 ا‎ 
۲۲ ) 1۲ + ۲ (۶۲ ۳ ([-)+7 -$y/Y9) 


r A 


) 1-1( 0 5 2 ری‎ | u ን. 022 9592672 
(1+ፀ )(1-5,/ሃሬ፥5(5.,/“7ሃሬ)፤፤ 1-5ህ /ሃ65፥፡(፤ 5» “ሃሬ? ' 
and we have used a/wáx*“za/wx=1/kx. lt is not clear from this 
relation how COP will vary with Prandtl number. However, if 
we define a set of variables and perform a numerical 
analysis we can see the desired effect. We will define a 
quantity, b4=( 5) 5.) constant; as the geometric mean of the 
penetration depths, such that ET and EY ELE Y IO. 
This roughly corresponds to our experiment. Setting 
reasonable values for the four parameters in our relation: 
وم عينم‎ Oe 5 Y 5 and Y=1.667, we can see how COP 
varies with Prandtl number. Figure [I-2 is a plot showing 
how heat (H5), WOrk (2) and COP vary with Prandtl number, 
using Equations II-46, II-47 and II-48. It is clear that as 
Prandtl number is reduced COP increases. It should be noted 
that as. Prandtl number becomes very small we actually see 
COP decrease with reduced Prandtl number, aS we estimated 
using the approximation with only the lowest order 
viscosity terms. Since Swift’s engine (1989) with CR 
sodium has a Prandtl number of approximately 0.004, his 


approximation is in fact valid for his system, although it 


is not for ours, as we discussed above. 
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b. BINARY GAS MIXTUURE THEO 


۳ sound Speed Theory 


One of the parameters of importance, when consi- 
dering refrigerator performance, is the-velocity of Som 
in the fluid being used. The purpose of this section is to 
show how this value is calculated for binary gas mixtures. 

In general, the sound speed "a" in a gas is given by 
the relation: a=(YR¿T/m) ”, where 52%/ና, is the 3۵ ۰ ۱ ۳ 


isobaric to isochoric specific heats, Ro is the universal 


gas constant, T is the absolute temperature and m is the 
molecular weight of the gas. This relation holds for a 


binary mixture of gases as well. However, Yandm may be 


- 


different for a mixture and must be calculated. Frome 
(1976) we get a relation for the specific heat capaci 


à binary mixture: (c_), =[2ኖ2ኛ. 8 ር ተ Mc ,.1/Dx, MOS ጽር 
P ij i ipi j j Pj i 1 e 
where X. is the mole fraction of the ps Spelexe imm. H iS NE 


, th 


molecular weight of the species and C pi is the isobemme 


te 
specific heat capacity of the i 5 Species. ۸ similar 


relation holds for c, as well. The equivalent molecular 
weight m, of the mixture is given by: ا‎ 
Combining the above relations we get 2 112 1 E 
ion for the speed of sound of a binary mixture, as follows: 
3 
a-(ERAT(x,M,c, ኀ48 ሸር C(x NM x mec 
0 1 1726.) 1. ከ) 1 ET i i Vi ከ j ۵ 
Since Y is approximately 5/3 for any ideal monatomic 


gas it is clear that any vario tao e Ca and C will be 


small. Therefore, the primary effect on sound speed, of 
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00 two Frases 15 that it 1S inversely proportional to 
E 252000516 TOOL “Of the total mass. At 300 degrees Kelvin 
the sound speeds of the mixtures used in our experiment are 


as follows: 


pure helium a*1008 m/s 
helium-argon (19% Ar) a*614 m/s 
helium-xenon (12% Xe) a*460 m/s 


2. Prandtl Number Theory 


a. General Discussion 
Before going into the details of how Prandtl 
number varies for binary gas mixtures, it is important to 
understand the general concept of Prandtl number, and to 


TUS መከር expected Variation based on this concept. 


randt l number is a dimensionless ratio that relates 
ከ ን በ፲[1፡ properties of fluids. In Fluid Mechanics it is 
given the symbol "Pr" and in Thermodynamics we use "vs", 


Prandtl number can be represented in many different ways. 
For calculational purposes we define Prandtl number as 
ን where N" is dynamic Viscosity, iS the isobaric 
Peeeit1e = near capacity and A is the thermal conductivity of 
the fluid. This relation is used because these parameters 
are well tabulated and more readily calculated than others 
Piet are available. FOr an understanding of the physics of 
Prandtl number and how it “relates to the refrigeration 
process we define Prandtl number as r=V/k, where v=VY/P is 


the kinematic 971506 517 and S عت‎ is the thermal 
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diffusivity. This definition is used because kinematic 
viscosity represents momentum diffusion through the fluid 
due to a gradient in the velocity, and thermal diffusivity 
represents the diffusion of thermal energy due to 
temperature gradients in the fluid. This says that Prandtl 
number is a ratio of momentum diffusion to thermal 


diffusion. A third way of describing Prandtl number is 


Cor 


i=(5,/5,)*, where v is the viscous penetration depth ana 
56 35 the thermal penetration depth. The capacity for 


diffusion of a given quantity is directly related to the 


magnitude of the associated penetration depth. Since 
A is a function of kinematic viscosity, and 
TESTA) is a function of thermal diffusivity, we will 


restrict the discussion to v and « effects. 

Hirschfelder 5 al., (1954, p.14) derive a 
relation, using simple kinetic theory, for the thermal 
conductivity of a monatomic gas, A=15R 9/48, where Ro 15 
the universal gas constant and WM is the molecular weight 
ofthe gas. From this relation and the knowledge that Co 
Cy =R and ۷ where R=E¿/M, it is easily shown that 
lS o ا‎ for a monatomic gas. It is also known 
that # =5/ሥ3 for any monatomic gas, which ር in በር. 
for all monatomic gases. Calculations using rigorous 
kinetic theory also show that Prandtl number is approxi- 


mately 2/3 for monatomic gases. This theory will be 


discussed in the next section. 
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When two monatomic gases are combined the 


molecular interactions produce changes in the physical 


properties, resulting in Prandtl number deviating from the 
2/3 value. For a mix of a heavy gas and a light gas, as we 
have used, the Prandtl number drops rapidly ۰ the 


ከዱ ርር Of heavy gas increases írom zero. This is due‏ ان 
primarily to a very rapid rise in the molecular density‏ 
=grams/mole) of the gas as the concentration of the‏ ( 
heavier gas increases from zero to some small percentage.‏ 
AS ah example, a mix of two percent xenon With 98% helium‏ 
produces a .64% increase in this molecular density. Since‏ 
the kinematic viscosity, V="/¢9, is inversely proportional‏ 
this molecwlar density, while dynamic ‘viscosity is‏ 
relatively invariant, there is a rapid drop in kinematic‏ 
viscosity and therefore momentum diffusion. Thus, we see‏ 
that the momentum diffusion process is dominated by the‏ 
presence of the heavier gas molecules. The thermal‏ 
diffusivity, AS Sn the ihe hand is inversely‏ 
to E Product fep Which is essentially‏ 29210221 05 ۳۳ 
constant With Changes IN Concentration, And directly‏ 
proportional to the thermal conductivity, thereby changing‏ 
due tom variations in thermal conductivity. Reif‏ 1217 01 
represents thermal conductivity by A=c(kT/m) “, where‏ )1965( 
CNMENS CG He capacity per molecule- (which is essentially a‏ 
eonstant), k ais the Boltzmann constant, T is the absolute‏ 


temperature and m is the molecular weight. It is clear from 
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this relation that the thermal conductivity, and therefore 
thermal diffusivity and thermal energy diffusion, will be 
much greater for a lighter gas, and their values will be 
dominated by the concentration of the lighter gas in the 
mixture. Combining these tuo effects it is clear that both 
momentum and thermal energy diffusion rates decrease with 
increasing concentration of the heavier gas. However, the 
rate of change in momentum diffusion is much more dramatic 
at the lower concentrations and drops off rapidly as the 
concentration of the heavier gas increases, while the rate 
of change of thermal energy diffusion is much more gradual 
and continuous over the concentration range. This results 
in the ratio of the two, which is effectively the Prandtl 
number, dropping rapidly from the 2/3 value for the pure 
light gas as the heavier gas is initially introduced, 
reaching a minimum as the momentum diffusion steadies out 
and rising back to the 2/3 value as energy diffusion 
continues to decrease while the heavier gas concentration 
approaches 100%. 
b. Kinetic Theory and Quantitative Analysis 

The purpose of this section is to discuss the 
theory, Uith its assumptions and limitations, and the 
equations used to calculate the Prandtl number of a binary 
gas mixture. Much of the information presented hes 
reproduced in Appendix A in order to allow the removal and 


use of Appendix A without reference back to the main body 
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01 1-5 > ፡ኖሞ፡፡ች ር  ” ” ”ግ ም ም 1 9۱16 9 ا‎ al. (1954) introduce a basic 
kinetic theory in Chapter one, develop a rigorous kinetic 
ቤቢ in Chapter seven and present equations for the 
calculation of thermal conductivity and dynamic viscosity 
roi nary gas Mixtures in Chapter eight of their text. 
Using these relatwemsmewith anothers for the isobaric 


Specific heat capacity of the mixture we can determine 


Prandtl number, ከኸ. ከወ ant reproduce "tre work of 
Hirschfelder, et al. (1954) we will present a discussion of 
the methods they used in deriving their equations, 


including the assumptions and limitations that apply. Then 
we will present all of the equations necessary for the 
Aita tion of Prandtl number for a binary gas mixture, and 
a sample cuve of Prandtl number as a function of 
concentration that was created using the output data of the 
program presented in Appendix A, which uses these 
equations. 

The development of the kinetic thoery of gases 
is based on knowledge of the distribution function that 
represents the number of molecules of a specific species, 
MR uniteelLumefr- element, about a point in space, with 
0 |2521 iss in a unit=range about a specified velocity, at a 
given Ms tant اف‎ tine. This “distribution function is 
represented by f, (FV tt). Hirschfelder, et al. derive the 
۱ ۲ 2131 ۳ —equatien. in terms Of this distribution function, 


which is the basisefor discussion.of.transport properties 
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of gases. The Boltzmann equation is given as: 


M ۵۶۱۱ 1/ at, 
UV + — | re ( e DE bdbd ¿dv, (112% 
ده ۱ هه‎ n, VU, ||. D J 


where the superscript "1" denotes a post-collision value, 
5፡1. is the 6-01 131510۳ initial relative velocity, bis 
the impact parameter, m. is the mass of the nn Species, X 


is an external force on the molecule, V is the velocity of 
the molecule and ና is a rotational impact parameter. The 
hydrodynamic equations of change (conservation of energy, 
conservation of momentum, etc.) are derived from the 
Boltzmann equation. 

Hirschfelder, ፎቴ al., use a perturbation method 
of successive approximations, which they attribute to 
Enskog, to obtain solutions to the Boltzmann equation. 
Using this method they derive a form for the perturbation 
function #, that depends on space and time through the 
quantities, number density, mass average velocity and 
temperature, and their spatial derivatives: 

(j) 


/ðr)+nE (C a). (11-50) 
j i J 


(1) 
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The coefficients A B and SC 


1" ; 1 each represent rather 


involved integral equations uhich must be solved. They are 
Solved by a series approximation method, using the infinite 
sonine polynomial series. The 16 polynomials are 
defined as: 

(m) 
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Linear combinations of a finite number of these polynomials 


a0 


are shown to be the solutions to the integral equations for 
A, b, and ope However, the transport coefficients 
(which is what we are after) are expressed in terms of the 
50112 polynomial expansion coefficients. These Sonine 
polynomial expansion coefficients involve complicated 
combinations of integrals, called "bracket" integrals. 


Chapman and Cowling (1939) express these bracket integrals 


as linear combinations of a set of integrals: 
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(11-52), where is the reduced mass of the colliding 


Hay 
molecules i` and j, e is the reduced initial relative 
speed of the colliding molecules, b is the impact parameter 
and X is the angle of deflection of the molecules. 
Hirschfelder, et al., have tabulated values of these 2 
integrals for frequently encountered bracket integrals in 
Table I-M of their text. Hirschfelder, et al., go on to 
define the transport coefficients in terms of these qii s) 
integrals. It is these relations that ue have used in our 
determination of Prandtl numbers for binary gas mixtures. 
Before presenting these relations it is impor- 
tant to understand the assumptions used in this derivation 
and the limitations they impose, if any. Since the usual 
definitions of transport coefficients only apply under 


conditions only slightly displaced from equilibrium, our 


results are limited by this restriction. 
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Only binary collisions are considered, which 
limits us to densities low enough that three-body 
collisions are unimportant. The theory predicts that 
viscosity 15 independent Of pressure at constant 
temperature. Hirschfelder, et al., state that experiments 
on nitrogen have shown that only a 4% deviation in 
viscosity occurs as pressure is varied from 1 to 60 
atmospheres. Our operating condition of 10 atmospheres is 
well within this range and should therefore closely conform 
to the theory in this respect. 

The theory limits itself to Classical Mechanics 
and thereby precludes low temperature situations where 
quantum effects become significant. These quantum effects 
are less than 1% for helium above 200 Kelvin (and even 
smaller for heavier isotopes). The coldest temperature we 
have achieved to date has been 193 Kelvin, under a no load 
condition, so our normal operations are all above 200 
Kelvin and we are not concerned uith quantum effects. 

The Chapman-Enskog (1939) theory uses a first- 
order approximation in their perturbation method solution 
to the Boltzmann equation. This is only valid when the 
gradients of the physical quantities (density, velocity, 
temperature, etc.) are small. This means the changes in 
these properties over one mean free path must be small with 
respect to unity. Based on the mean free paths for gases at 


pressures above one atmosphere, this condition is satisfied 
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Seo under eendi tions otr extreme gradients, such as those 
present in E waves, Since the observed pressure 
waveforms in our resonators are sinusoidal, we do not have 
any Shock waves and we are not restricted by this 
assumption. 

uc ፣፡፣ ther assumes that the dimensions 
of the containing vessel and any obstacles within it are 
large compared to the mean free path, which provides no 
limitations for our gases and pressures. Although the 
Chapman-Enskog theory applies strictly to monatomic gases 
the results . have been shown to be good even for polyatomic 
gases, provided that the molecules are not too non- 
spherical. Finally, the development of the theory results 
TS Set of equations for transport coefficients which rely 


upon ۲۱ لم‎ 19۵ discussed رن‎ መሌ 


integrals" These 
115551215 are dependent upon the potential function that 
ሱፕ” መ በቲ= the nolecular interactions. Throughout the 
derivation by Hirschfelder, et al., and in all of the 
ከ OMS “secar fom our caleulations, the Lennard-Jones (6- 
12) potential is assumed. 

NN تا‎ 9 ۱ e aed the theory and its 
limitations we will present the equations necessary for the 
9۱ ۳۸ On «of dynamic Viscosity, thermal conductivity and 
ዝ ፣ ፥፣ውሽኝጅክጅጂ፡ያዯሪሪ፡ቹ፡፡ heat capacity of a” binary mixture of 
In uM cases. Fron these, "as We håve already discussed, 


MEL set Pranatl number, where A 
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The dynamic viscosity for a binary gas mixture 
is given by: 
1+4, 


aix Ji s ; 
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| [ 


and 
3 M MPA o Sor ERE M M 
2 NL a SOR بجعت ,یف )نی‎ el ። 
3 As 4MM, (NJ, nad, ዘ) 
In these equations Xy and Xo represent the mole fractions 


of species 1 and 2, M and M, represent the molecular 


1 2 
ቀ 
weights of species 1 and 2 and ል; 7 11) 1)ቀ 


the superscript $* refers to a value that takes into account 


, Where 


the deviation of a particular molecular model (in our case 
the Lennard-Jones (6-12) potential) from the idealized 
rigid sphere model of molecular interactions. The values 
for C74, and a3, are given by the equation for a pure 
substance: 


(1,x10/-266. 93 (T) Pt e a6 0 £1 


UR 

where n is dynamic viscosity in gm/cm-sec, T is temperature 
VES ፐ”=ዜፕፐ/ ‹ is reduced temperature, HM is molecular 
weight, ሆ iS collision diameter in angstroms and :/k is 


the potential parameter in ?K. The remaining term ۲۷ و ولو‎ 


represents an artificial quantity for a hypothetical pure 
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Substance with molecules of molecular ueight 24,4,/(4,+B5), 
that interact ‘according to a potential curve specified by 
interaction parameters 712 and 12 These terms are 
described in terms of the Lennard-Jones (6-12) potential, 
where ያ is the distance of closest approach (in angstroms) 
of two molecules which collide with zero initial relative 
kinetic energy, and ፥ is the maximum energy of attraction 
of the two molecules. (12? 1 is given by: 

۲1 مر‎ 1۵12266۰936211 M, T/ (M 48,0 3 "7C e S20 7 P f op, 


where the subscript 12 on individual terms refers to a 


combined value for Species 1 and 2. The values for 712 and 


C12 are given by the empirical combining laus )2۶و ؟‎ ۲ ( 
B. 2 | X ; 
and C157 Cf4 65) i respectively, and Ti2 is given by 
m ሩ 
፲12 -kT/ 5,5. 


The thermal conductivity for a binary mixture 


is given by: 
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The definitions for all terms are similar to those given 
for the dynamic viscosity relations. 

The remaining quantity necessary for the calcu- 
lation of Prandtl number is the isobaric specific heat 


capacity of the binary mixture. Huang (1976) gives a 


: = 2 3 £ 
relation for this quantity Comix mc pi /m- ES x,H, 
Where X, is the mole fraction of species i, M, is the 
molecular weight of species i and C pi is the isobaric 


Specific heat capacity of species i. The required units of 
this quantity, for compatibility With the previous 
equations are cal/gm-?*K. 

All of these equations were formed into the 


Fortran computer program presented in Appendix A. Outputs 
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have been created (from this program) for temperatures of 


2004© and 300*K, and for mixtures of helium-xenon and 
helium-argon, for all concentrations from 100% helium to 0% 
ከእጅ ከ AN composite plot Showing all of these cases is 


Pad که‎ ۱۱۶ LI-S. It is clear from these curves that 
ን “ስ ተ 111 2۳ number variation with concentration is just as 
predicted by the qualitative analysis, given in the first 
ECE Of this section. 
3. Selection of a Gas Mixture 

Now that we have shown how the sound speeds and 
Prandtl numbers of binary gas mixtures vary with concentra- 
tion, it is necessary to consider how each of these 
parameters affect performance. We know from the 
thermoacoustic theory, presented in section A of this 
Chapter, that power density will decrease with reduced 
sound speed (due to the increased mass), and that COP will 


increase With reduced Prandtl number. However, as Prandtl 


number decreases, due to the addition of a heavy gas, the 
ewe ۱۱2 92 the gas ۱ ۱16 2 5 6 5, and sound speed is 
reduced. It is therefore necessary that the mixture have as 


lou a Prandtl number as possible while the mass increase i 


LA 


kept to a minimum (and therefore the reduction in sound 
لد دد زد‎ minimized). Locking at Figure II-3 it is clear 
that we need a mixture that iS at a point on the very 
Stecply SDN section of the curve, as the concentration 


begins to deviate from pure helium. Our selection of 12.5% 


SN 


He-Ar 300 
He-Ar 200 
He-Xe 300 
He-Xe 200 





Helium Fraction 


Figure [1-3. Prandtl Number vs. Helium Fraction 
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Prandtl Number 


Bund CES 31 201 if our helium-xenon and helium-argon 
mixes, respectively, was made based on this requirement. 
With the xenon we achieved approximately 90% of the maximum 
possible Prandtl number reduction, while getting only 40% 
Ie possible Mass increase. With argon ue once again 
achieved about 90% of the maximum possible Prandtl number 
reduction but got approximately 60% of the possible mass 
increase. The results of these selections will be discussed 


NAO napter V. 


C. STACK DESIGN THEORY 

MINAS Of this Section iS tO explain the diff- 
erences between the tuo stacks built for this project and 
mat used if previous experiments, including the basis for 
any and all modifications. To provide an orderly flow of 
information, the differences between the stacks will be 
fully described and then an explanation of each of the new 
features will be given. Details of stack construction are 
provided in Chapter III. 

 “ ل‎ |1701 three Slacks: in this section: the 
tema Stack, built by Dr. Tom Hofler, will be referred 
eee ord Stack, The first Stack built for this 
experiment will be referred to as the "new" stack, and the 
second stack built for this experiment will be referred to 
PEN un Stack. The actual designs for all of the 


Stacks are based on the work and ideas of Dr. Hofler. All 


E» 


of the stacks have the same outer diameter of 1.506 inches 
and the same size center post with a 0.25 inch diameter. 
Additionally all three have the same plate thickness of 
0.003 inches, and were built using a loom, discussed in 
Chapter III, with a line spacing of @.2 inches. The primary 
differences between the stacks are in the areas of average 
plate spacing, overall stack length, and final arrangement 
of the fishing line used to provide the plate spacing. 

The old stack has an overall length of 3.089 inches with 
an average plate spacing of 0.0140 inches, with no 
modifications to the fishing line arrangement. The new 
Stack has an overall length of 3.9066 inches with an average 
plate spacing of 0.0133 ا‎ Additionally, it has a 
Special pattern cut into the fishing line, as shown in 
Figure II-4. The dual stack has an overall length of 3.068 
inches with an average plate spacing of 0.016 inches at the 
hot end and 0.035 inches at the cold end. There sS 
pattern cut into the fishing line of the dual 
however, there is a gap of approximately 0.75 inches 
between the rows of thin lines and thick lines, as shown in 
Figure 11-5. Additionally, on the layer of plastic ss 
both sizes of fishing line, every other piece of thE በዘር 
fishing line has been removed, starting at a point 18 
inches from the end of the plastic sheet that is at the 
center of the roll. The only remaining difference ۶ 0 


the stacks is that the old and new stacks have center posts 
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Ii COTS 6110-0 type rod and the dual stack has a nylon 
ED cad lhrs last difference is strictly due to 
availability of material at the time of construction and is 
not considered to be of any Significance in relation to 
refrigerator or stack performance. 


Each of the modifications, relative to the old stack, 


will now be explained, beginning with those for the new 
Semen. The changes in this stack are primarily intuitive in 
nature, based on experiments by Hofler and Wheatley, and 


therefore do not lend themselves to quantitative analysis. 
The change .in overall length was selected to provide a 
Space between the end of the stack and primarily the cold 
heat exchanger. This space is desired to preclude any flow 
mere trots = that COULd CCCUr With the Stack plates resting 
2551151 the plates of the™heat exchanger surface. A piece 
UT fishing line (approximately 0.020 inches in diameter) 
was glued to the end of the stack, to fill the void left by 


the reduced length, thereby preventing the stack from 


moving around in the resonator relative to the heat 
0 252112515. The pattern of Figure II-4 uas selected for tuo 
MIT suce the amount of heat transport is related to 


the total surface area of the plates, the removal of some 
Of the fishing line would increase this surface area and 
IM case” the heat transport. Also, since it is 
راد‎ 0012-1  inpossibie to have precisely uniform spacing of 


۱۳ ۱ 115112 line and plates, there are slight differences 
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in the longitudinal flow channels in the Stack. DAN 
turn can produce flow variations and pressure differen mis 
from one channel to the next. Due to there being no 
continuous pieces Of fishing line along the full lens OI 
the stack, With this pattern, the channels are all 
interconnected and the pressure and flow differences, if 
any, will be minimized. The overall result of this pattern 
was a reduction of approximately 50% in the amount of 
fishing line present in the stack. This reduction should 
also have the effect of reducing the overall obstruction sS 
flow in the resonator. lt should also be noted thal sS 
pattern may have had à significantoeffect on the deum! 
crush experienced during the rolling process, resulting in 
the smallest average plate spacing of any of the stacks. 
This plate spacing is very important in the thermoacoustic 
performance of the refrigerator. 

The changes in the dual stack are somewhat nore ON 
and theoretical in nature. Within the accuaracy of 
construction methods, the overall length of the dual stack 
is the same as that of the new stack, for reasons explained 
above. The removal of every other piece of the thick 
fishing line, after the fir l S man እጅ ከ is based on 
reducing the overall amount of fishing line. in the stack; 
as discussed for the new stack. All lines were left oum" 
for the first 18 inches to provide a more uniform pihake 


spacing toward the center of the roll, where the diameter 
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of the roll is not very large with respect to the fishing 
line spacing. The 0.75 inch gap between the rows of thin 
and thick strings was used for two reasons. A certain 
minimal gap iS necessary to prevent the inner plastic 
layer, on the hot end, from butting up against the thick 
Strings when the ends of the stack are pushed together. The 
extra distance, out to 0.73 inches, was used toa give 
approximately the same total length of string as we had in 
the new stack. 

Understanding the reasons for the dual nature of the 
Stack requires some knowledge of the basic characteristics 
of large temperature span refrigerators, as well as the 
basic heat transfer process that occurs in the stack. The 
heat flow at the two ends of the stack is not the same, 
generally being significantly higher at the hot end, where 
heat is rejected. This difference is functionally dependent 
on efficiency and increases with reduced efficiency. To 
properly explain this effect some terms must be introduced. 
Although the following relations are simplified they are 
adequate for the explanation of this effect. The heat 
transfer capacity q is given by: 

ባ(፳%)=9(፳)/(የቦ(፳)-1)=-5,5(፳)9(5), (11-52) 
and the local heat flou Q is given by: 
Q(x)z- lá, p(xJ)u(x)(Cr(ix)-1), au 539 
where r(x)=0T/0T vT is the maximum temperature 


crit Crit 


gradient in the x direction, Tis the surface area per unit 
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length of the stack, p(x) is the acoustic pressure 
distribution, u(x) is the acoustic velocity distribution, 
5.4 is the thermal penetration depth (which according to 
Swift (1989), is the distance that heat can diffuse through 
the fluid during the time 1/w), and x is the longitudinal 
distance along the stack. In a standard stack u(x) and f(x)- 
1 are the only parameters that vary significantly in the 
longitudinal direction. The quantity nr(x)-1 represents the 
temperature gradient perpendicular to the x direction (i.e. 
the transverse direction). This quantity will vary greatly, 
depending on heat load, and must be large in order to 
transfer large amounts of heat. In a large temperature span 
refrigerator የ(፳)--1 m very small at the cold end and 
large at the hot end. The Q(x) distribution is dominated by 
|roo-i[. Conversely, u(x) is generally small at the hot 
end of the stack and large at the cold end. It is u(x) that 
dominates the q(x) distribution. These effects give a low 
heat flow at the cold end, where the heat transfer capacity 
15 highest, and a high heat flow at the hot end, where the 
heat transfer capacity is lowest. 

The object of the dual stack design is to reduce this 
conflict by changing the distribution of q(x), making it 
more uniform and similar to that of Q(x). Since nr(x)-1 must 
be geater than zero to transfer heat, and is very dependent 
on heat load, it was determined that q(x) was the quantity 


that should be modified. To accomplish this, the acoustic 
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velocity distibution had to be countered by some other 
effect. The only parameter affected enough by stack design 
to accomplish this is V. By varying the amount of surface 
area along the stack length the q(x) distribution can be 
altered. By having half the number of plates at the cold 
end, we cut the surface area by a factor of two, 
essentially countering the u(x) distribution, and making 
the capacity for heat transfer correspond more closely to 
the actual heat flow conditions. Ideally, the surface area 
Should vary continuously in the longitudinal direction to 
precisely counter the velocity distribution and produce an 
exactly uniform q(x) distribution, thereby giving optimum 
results. This situation would be much more difficult to 
construct and methods that have been discussed to date have 
other drawbacks that make them less acceptable than the 
current design. Although this is the major reason for an 
expected increase in the COP with the dual stack design 
there are two other minor effects which should also help. 
An area not covered by the simplified equations above is 
that of viscous friction as the gas moves along the plates. 
The energy loss due to this friction is proportional to 
both the square of the acoustic velocity and the total 
Surface area exposed to the gas. Since the surface area has 
been cut in half at the cold end, where the velocity is 
highest, this source of energy loss should be significantly 


reduced, thereby increasing the COP. The second area of 
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interest relates to the uniformity ot the q(x) 
distribution. Since q(x) depends on u(x), anything that 
Makes u(x) more uniform will also make q(x) more uniforn, 
Which in turn tends to improve the COP. The velocity 
distribution will depend somewhat on the degree of blockage 
in the flow path. The plates have an "effective" size that 
is larger than their thickness, due to viscous effects. By 
reducing the number of plates the degree of blockage is 
reduced and the velocity distibution is more uniform, 


thereby tending to increase the COP. 
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III. STACK CONSTRUCTION 


The thermoacoustic stack, necessary for the transport of 
heat as discussed in chapter II, consists of a continous 
roll of very thin plastic, the layers of which are 
separated by spacers attached to the plastic surface. This 
geometry was selected instead of a parallel plate 
arrangement based on both ease of construction and 
structural rigidity of the finished product. When the 
dimensions of a stack are considered, plate spacing of 
0.015 inches (or 0.033 inches depending on which stack), 
plate thickness of 0.003 inches, overall diameter of 1.506 
inches, and overall length of about 3.066 inches, it is 
easily seen that the parallel plate geometry would be much 
nore difficult to construct and extremely fragile. The 
plastic used in this experiment was a polyester Mylar TM 
` type. There was no particular reason for this selection 
other than ease of handling. The plate spacers consist of 
strips of high grade monofilament nylon fishing line which 
are affixed to the plastic surface using a spray adhesive. 

Prior to construction, decisions must be made concerning 
overall length of the stack, which determines the width of 
the plastic sheet, required diameter, which determines the 
overall length of the plastic sheet and, approximate plate 


spacing, Which determines the size of spacers to be 
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attached to the plastic surface. For this project two 
Stacks were built. The first had an overall length of 3.066 
inches, with a diameter of 1.506 inches and 0.015 inch 
spacers. The second was a dual stack with the same 
diameter, an overall length of 3.9068 inches and spacers of 
0.015 inches, with two layers of plastic along half of its 
length, and spacers of 0.033 inches along the other half. 
This provided a variable plate spacing for reasons 
discussed in chapter II. The overall length of plastic 
sheet needed to make a stack can be determined by 
mathematical analysis of the spiral wound roll once the 
dimensions are specified. The resulting equation for length 
is L=(t/d)(r,*-ro?), where dis the sum of the spacer size 
and plate thickness, r, is the radius of the stack (and 
also the radius of the resonator bore), and ro iS the 
radius of the center post, on uhich the plastic is rolled. 
The center posts that have been used are 0.25 inch diameter 
and are either cloth phenolic or a nylon type plastic rods. 
Usually a safety margin of twelve to eighteen inches is 
added to the length calculated by the above equation. This 
is due to the fact that a precise fit is desired in the 
resonator bore and it is easier to trim small pieces from 
an oversized roll than to add pieces to one that is 
undersized. The width of the plastic sheet is cut very 


precisely using a sharp paper cutter, with a guide that has 


been carefully positioned to the desired width. This method 
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has provided results with an accuracy of 0.002 inches. The 
overall length of plastic sheet used for the two stacks in 
this experiment vas 120 inches for the first stack and tuo 
pieces of 60 inches each for the dual stack. 

Once the plastic sheet has been cut to size and fishing 
line cof the required size has been acquired, construction 
can begin. The spray adhesive used in this experiment uas 
3M brand Super 77 spray adhesive. Since attaching 
individual pieces of fishing line would be very time 
consuming and result in a less than uniform spacing an 
alignment device is used to assist in this process. An 
aluminum loom has been constructed to provide uniform 
parallel spacer positioning on the plastic surface. A long 
piece of fishing line iS wound around the hooks of the 
loom, positioned in grooves to provide parallellism, and 
fixed at both ends with screws to hold it in place. The 
loom has an open front and back to provide access to both 
Sides of the fishing line. Using a spray shield to protect 
the loom, spray adhesive is applied to the front side of 
the fishing line. The loom is then turned over and aligned 
over the plastic to provide parallel, uniform positioning 
of the strips of fishing line. Once positioned the loom is 
set down and a smooth aiuminum block is used to press the 
fishing line onto the surface. The working surface and 
press used in this experiment were smoothly polished 


aluminum blocks, which vere necessary to provide a 
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relatively uniform compression of the fishing line when 
pressing. A teflon sheet was placed under the plastic and 
over the working surface to prevent the strings, which 
extend beyond the plastic, from sticking to the working 
surface. A sketch of the loom, with fishing line attached 
is shown positioned over the plastic sheet in Figure III-1. 
The distance between strips of fishing line on the loom 
used in this experiment was 0.20 inches. This spacing was 
selected based on experiments conducted by Dr. Tom Hofler 
at the Naval Postgraduate School. The loom covers a length 
of six inches and therefore requires twenty repetitions of 
this procedure to complete a stack. The dual stack required 
more repititions due to there being two different string 
sizes on one plastic sheeet and a second layer of plastic 
of half the width, as shown in Figure III-2. Once the 
fishing line is pressed down and allowed a minimal drying 
time of one to two minutes a Scalpel or razor blade is used 
to cut the fishing line off at the edge of the plastic. It 
Should be noted that if the line is wound on the loom with 
too much tension the plastic will tend to curl when the 
line is cut. This is obviously an undesireable situation so 
care must be taken when winding the line on the loom to 
provide just enough tension to maintain the strips in the 
grooves. After all of the fishing line is attached to the 
plastic some trimming must be done prior to rolling. The 


fishing line is carefully cut back approximately one 


ue. 


Figure III-1. Stack Construction Illustration 
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Sixteenth of an inch from each edge, to prevent potential 
flow problems at the heat exchanger interfaces. 
Additionally, depending on the stack design, a pattern is 
cut in the strings or some strings are removed, for reasons 
explained in the stack theory of chapter II. Great care 
Should be exercised in the trimming of any patterns to 
avoid making any large cuts in the plastic surface as this 
could produce Significant flow variations during operation 
of the unit. 

After completing this process the stack is ready to 
roll. The center post is carefully taped, using scotch 
tape, to the end of the plastic sheet. The plastic is 
carefully rolled up, using just enough tension to maintain 
the Shape and provide uniform layer spacing, without 
getting too much crush, which would result ina layer 
separation much smaller than the string diameter. This is 
in general the most difficult part of the stack 
construction process. The uniform layer spacing is checked 
visually by holding the stack up to a light or using some 
sort of magnifying device. The spacing can also be checked 
quantitatively based on the known innner post diameter, 
the measured outer diameter and a count of the number of 
layers along several radii. This will provide an average 
layer spacing which will in general be somewhat smaller 
than the sum of the fishing line diameter and the plastic 


thickness. This is due to a certain degree of crush or 


5 


compression during the pressing and rolling processes. 
After rolling, even With a great deal of care, there will 
be some degree of spreading, or misalignment, at the ends 
of the cylinder. Using two smooth surfaces the ends of the 
stack can be pressed and the layers will slide into 
position to give flat ends and an overall length as precise 
as the width of the plastic sheet. 

Drawings of the two ends of the dual stack are shown in 
Figures III-3 and III-4. The other stack has the same 
appearance as the dense or hot end of Figure III-3. 
Polaroid photographs of the actual dual stack, showing a 
radius from the center post out, at a 6.4X magnification, 
are provided as Figures III-5 and III-6. The finished roll 
is then trimmed down in diameter, by cutting small pieces 
Off the length of the plastic strip, to fit snugly in the 
resonator bore. The final product is taped along the sean, 
With scotch tape, to hold the roll in the cylindrical 
shape. 

Although some of the differences between the two stacks 
that were built for this experiment are mentioned in here, 
the details of and reasons for the differences are 


presented in the stack theory of chapter II. 
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Figure III-S. Photograph of Dual Stack, Cold 
End, at 6.4x Magnification 
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Figure III-6. Photograph of Dual Stack, 
at 6.4x Magnification 
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IV. EXPERIMENTAL APPARATUS AND PROCEDURES 


In this Chapter we will present a description of our 
experimental apparatus, including the electronics used, and 
a discussion of our experimental procedure. We have used 
the prototype thermoacoustic refrigerator, and associated 
electronics package, that was built by Hofler (1986, 13887. 
We have paraphrased sections of his” discussion በህ በ 
system as necessary to provide a basic understanding of the 


apparatus and its operation. 


A. 12150111 1 15 5 APPARATUS 

in order to properly describe the apparatus 8 
necessary to break it down into two major areas; that of 
the high intensity driver, With its instrumentation for 
acoustic power measurements and that Oof the refrigeri mena 


components (resonator, stack, heat exchanger, etc.). 


1. The Driver and its Instrumentation 
e driver consists of a modified commercial 
loudspeaker, shown in Figure IV-1. The modifications 


consist of removal of the original ۳128 ۲ 16 ۳0۱ 5 ل‎ ህ፣) 
fabric dome and attachment of a thin-wall aluminum cone- 
shaped piston and a surround. Figure [V-1 shows the unit 
With a voice coil and former 1, iron pole pieces 2, magnet 


3, aluminum piston 4, surround 5, aluminum pressure vessel 
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Figure IV-1. Schematic of Driver Apparatus Shouing 
Acoustic Motor, Transducers, and a Portion of a 
Resonator. Numbered Parts are Identified in the Text. 
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6 and 7, aluminum plate 8, cooling water tubes 9, annular 
plate (where pressure transducer is mounted) 10, dynamic 
pressure transducer 11 and a post for the mounting of an 
accelerometer 12. The unit has a two inch diameter voice 
coil and an effective piston driving area of 0.93 square 
inches. Additionally, it has a high power capability Ts 
is made possible by suspending ferrofluid in the magnet 
gap, to conduct heat from the voice coil. A capillary 
filled hole is located adjacent to the surround to break 
the gas seal provided by the surround. The flow impedance 
of the capillary was chosen to be as low as possible 
Without affecting the acoustic power measurements. The 
assembly is housed in the aluminum pressure vessel 6 and 7, 
and the resonator section is bolted to 6. The heat removed 
by the water cooling tubes ሠ, is conducted through the 
aluminum housing 6 and plate 8. 

The dynamic pressure transducer 11, is physically 
attached to the structure of the driver and can therefore 
be calibrated along with other aspects of the driver, 
independent of the resonator. The transducer 1 ۳ ۵ ደጀ ከርከር 
designed and built by Dr. Hofler. It consists of 3 2111 
Cut quartz crystal disk With evaporated electrodes and a 
miniature FET impedance converter circuit that buffers the 
Output signal. The dynamic pressure transducer provides a 
low noise, low distortion, phase accurate signal Which is 


calibrated to determine the volume velocity in a small 
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SION "Gavity that Cah be attached to the driver. 
Ns, NENNEN iON Withe=tme integrated accelerometer 


EDI (Which is then a velocity signal) and the phase 


difference between the two, provides for an accurate 
measurement OT the acoustic power provided to the 
03 0301. Hofler's calibration of the pressure transducer 


۲۳ 20000051151160 by comparison to a "standard" transducer 
that had previously been calibrated to an accuaracy of 11%, 
۱3 ۷11 ۲۷ et. al. (1982), resulting in a sensitivity of 
۱۳۳ 9 (11%) volts/bar. The accuracy of the acoustic power 
measurement is 13%, as shown by Hofler (1988). 

6). ..!ና ርስዕ ርጋዕስስም ክስር discussion of Figure IV- 
1, the accelerometer is mounted on the post 12, which is 
attached to the back side of the piston 4. The post is 
necessary to thermally isolate the accelerometer from the 
Piston (which becomes warm during operation) and to 
position the accelerometer where its cable has clearance 
from the piston wall. The accelerometer is a shear 
piezoceramic type. The piston’s velocity signal (which is 
the integrated output of the accelerometer) is calibrated 
Ea tie  sayilemic pressure transducer signal using a small 
nommer ال اک‎ ACOUSTIC impedance, to get an accurate 
value for volume velocity. This calibration may change as 
the properties of the suround change with use and must 
therefore be repeated at periodic intervals, to maintain 


the accuracy of results. Repeated calibrations, over a two 
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year period, have shown that the sensitivity of the 
pressure transducer has been stable within about ð. 17; 
while the volume velocity calibration value has shown an 
overall decrease of about 7%. It should be noted that the 
volume velocity value dopped just over 4% during the first 
four months of driver use, and dropped only 3% nore jeer 
the next 18 months. 
2. Refrigeration Components 

The main components of the refrigeration Section are 
Shown in Figure IV-2. The resonator (18, 20 and 21 MS 
figure) is constructed of several parts. Where tempera es 
ungformity is desired heavy-gauge copper is used: lan 
Section 18, containing the stack 17, there isa large 
temperature gradient and an epoxy-fiberglass substance is 
used to provide high strength with very low thermal 
conductance. This epoxy-fiberglass section, by itself would 
be very permeable to gas diffusion if left alone. To 
correct this situation, Dr. Hofler evaporated a thin film 
of indium over the outside of the fiberglass sectio m 
then wrapped more fiberglass över the meta E TE 

The primary refrigeration component is the stack ۲ 
in the Pigurer 1132) whieh is 56 of a long sher TEES 
plastic that is spiral wound on a- 0.25 inch dial sS 
plastic or cloth phenolic rod 16. The layers of plastic are 
separated by monofilament fishing line spacers. The detaiTs 


of stack construction uere presented in Chapter III. 


84 





21 


Figure IV-2. The Primary Refrigerator Components, 
Shown Attached to the Driver Assembly. The Numbered 
Components are Discussed in the Text. 
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The other components (eor inportance in tie 
refrigeration section are the hot and cold heat exchangers 
15 and 22 in Figuure I[V-2, respectively. The hot heat 
exchanger provides for removal of heat from the unit, to 
prevent overall heating of the stack during operation. The 
cold heat exchanger provides good thermal contact with the 
applied heat load, from the heater coil 19. Dr. Hotler 
constructed the heat exchangers using electroplating and a 
chemically removable aluminum form, shown in Figure 1۷-3 
prior to removal of the aluminum. The resulting Monk 
exchangers had parallel copper strips with the planes of 
the strips parallel to the tube axis. The details of this 
construction are provided in the doctoral dissertation of 
Hofler (1986). 

The remainder of the physical apparatus consists of 
a vacuum system, that attaches below the driver housing, 
thereby enclosing the entire refrigeration section, and a 
section of  super-insulation wrapped around the resonator. 
All of this is to provide thermal insulation, them 
reducing nuisance heat loads. A layer of super-insulation 
consists of a layer of glass micro- fimer web followed by a 
layer of aluminized polyester film. The entire rese 8 
Section is Wrapped with 10 to 12 of these layers. The 
vacuum system uses an ion gauge to monitor the vacuum, 
which varies, between Sml@n*. and .Sxieue foun during 


Operation, depending on operating conditions and gas types. 
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The electronics portion of our experimental apparatus 
consists of two primary sections. These are the control") 
protective section and the data monitoring section. The 
control and protective systems will be discussed here and 
the data monitoring aspects will be covered in the next 
section, where ue discuss the experimental procedure. 

The control and protective electronics system provides 
feedback control of the drive frequency and amplitude, 
protection against burn-up of the driver voice coi tea 
signal conditioning to increase the precision of acoustic 
power measurements. The protective circuit is particularly 
important in our situation because the required operating 
conditions are very close to the maximum limits of the 
driver. 

Figure IV-4 is a block diagram of the control CAMAS 
for drive frequency and amplitude. This figure is 
reproduced, With minor equipment updates, from the doctoral 
dissertation of Tom Hofler (1986). The explanation 013 ۶ 
Operations för the most parti is quoted here for 
completeness: 


"The frequency control is based on a lock-in amples 
that is used as a phase null detector with an 9 


Output. The velocity signal is fed into the DM 
preamp, and the pressure signal is used as a reference 
Channel input. The level of this reference signalis 
adjusted to be 1 V rms, in order to minimize lock-in 


phase errors. The Y-channel output voltage of the 106۳ 1 
is proportional to the measured phase deviation and us 
fed into a simple integrating controller. The controller 
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Output is connected to the ۷۵1۵86۰-001۲ 9 ۱ ۱۶0 ۶ 1 ۶ ٩ 
(VCF) input of the signal generator, and the loop is 
closed by the driver-resonator system. The phase, 
including calibration corrections, is set on the lock 
so that the desired Y-channel output voltage is zero when 
the pressure and velocity are in phase. The contro 
has a set point voltage fixed at zero, so the frequency 
is controlled to maintain the resonance condition of the 
pressure and velocity being in phase. The pressure 
amplitude control group begins by adjusting the amplitude 
of the pressure Signal to a pre-defined level with an 


adjustable gain stage. This ensures that the loop gain 
and stability margin are constant for different drive 
levels. After the pressure signal is adjusted mim 


amplitude, it is rectified and filtered, and this DC 
signal is used by another integrating controller ss 
produce a control voltage for an adjustable gain control 


Circuit (AGC). The output of the signal generator is then 
passed through the AGC and into the power ampliar 
input. The loop is again closed by the acoustic SyS Tens 


With its pressure amplitude thereby being regulated. 
Integrating the -accelerometer signal 1 57 ۰ 
since it is the velocity that is the important quom ا‎ 
in à power measurement...." 
To finish the explanation we will need to discuss the 
velocity signal measurement since the method used by Hofler 
was not used in our experiment. In our experiment we 
measured the velocity amplitude using the lock-in directly. 
It was determined that this would not introduce any errors 
into the power measurements, as previously thought. Unlike 
the previous apparatus, our lock-in amplifier uses true 
Sine wave conversion, resulting in complete rejection of 
the harmonic distortion in the accelerometer signal. Usine 
the measured dynamic pressure and velocity, along VIAS 
sensitivities of the transducer and accelerometer, and the 


volume velocity calibration factor we can calculan 


accurate acoustic power, as long as the frequency control 
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EM cL iS maintaining the phase difference between the two 
Signals at zero. 

۱۱۳ power Protection Circuit monitors the power amp- 
lifier current and uill trigger a crowbar circuit to blow a 
በ. either Of two setpoints are exceeded. lf a set peak 
current value is exceeded the fuse is blown, to protect 
۳۳۳۳۲ ۰ large fast transients. If a set rms current value 
is exceeded the fuse will also be blown, to protect against 
currents below the peak limit but high enough to burn out 


the voice coil. 


C. EXPERIMENTAL PROCEDURE 

"1 اد وا‎ We Will discuss Three basic areas of our 
experimental procedure. First is the setup of the apparatus 
for a specific data run. Then we will cover the procedure 
Iu starting up and operating the refrigerator. We will 
Ende With a discussion of the actual data that is 
۱۳۳ تا‎ in the operating section we will cover two 
different methods since the procedure was changed after the 
second data run, for reasons that will be explained later. 

before discussing the actual setup it should be noted 
A temperatures that we recorded were measured using 
commercial type-E thermocouple wire with an accuracy of 
t1.7*C or 11% of the reading, whichever is greater. All of 
the thermocouple wires are fed through a Hewlett Packard 


TANS data acquisition “Control unit to a computer, for 


Ed 


۵ 11 የ display on the computer monitor during 
Operation. 
1. Setup of Apparatus 

The setup of the system for a particular datos 
consists of filling the system with the desired gas, 
pressurizing the system up to the desired pressure of © 
bar and starting up the external vacuum systen. 

The resonator is filled through a connection adde 
top of the driver housing, however the controlling factor 
is the capillary near the driver surround that Was 
previously discussed. The pressurization rate iS monitored 
on a differential pressure gauge since the surround aie 
capable of sustaining no more than 0.25 bar of differential 
pressure without permanent damage. The system is filled 
using a process Of evacuation and purging that is Fees 
twice before the system is brought up to operating 
pressure. This procedure is designed to remove an EE EE 
Other than the desired gas (particularly air if the لمج‎ 
has been opened) prior to system operation. After the 
second evacuation the system is pressurized to 190 bar. The 
evacuation and purging process consumes the better 620 
à day due to the lou flou rates that are possible Within 
the alloued differential pressure range. 

To save time, the external vacuum system can be 
started in parallel with the above procedure. Since the 


diffusion pump of the vacuum system creates a great deal of 
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Heat, an external cooling supply is used to keep it ata 
reasonable temperature. Therefore, to start up the vacuun 
5 ۱ the cooling system is lit off, the main vacuum pump 
SWetion iS Closed, the vacuum pump is started and the pump 
Suction is slowly opened, tTO prevent rapid 
depressurization and damage to the previously discussed 
۳ ۱ 1 112 ۳1 8۰ Then the diffusion pump is turned on. 
Once the system is up to pressure and the vacuum is below 
1x103 torr it is ready to start up. 
ስስ ፡ ከ ብ Startup and Operation 
To start the refrigerator the tuo controllers 
Es ከጊ በር. and pressure amplitude) are turned off and the 
Signal generator is set up for a sine wave at the 
En -amate resonant frequency for the resonator and gas 
type being used. The drive amplitude is set to zero and the 
Bole! is turned on. The amplitude is brought up slowly 
until a drive current of about 9.5 amps is achieved. At 


ome Point the frequency is manually adjusted to resonance 


ከ፡ ፡ ጌስ quadrature output) and the frequency control 
Sere wit 15 turned on. Then the amplitude is slowly 
ከከ ስ ተት the desired value (Cif possible) and the 


reference channel is adjusted to precisely 1.00 volts. Then 
the 1م‎ ን ቲያልፕፕፕያሇያፕታፕቨ ውጤ ስጠጋ ከከህከርን Control 25220318 is turned on. 
 ኢሕኤሒይጅጂኪፏሞይጅ፣ፕሇ፡፡፡ ጠመ ጠጠ ፡ቶ፡. ር ሮም ን ርርዐ[ pressure amplitude and gas type it 
IU 6/,7፡.135515ሀ51ሮ ጊዕ reach the desired operating point 


Mimic startup. It is therefore necessary to select 
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some lower amplitude, until the engine cools doun and the 
desired amplitude can be reached. The limiting ۶ 
generally the piston’s displacement amplitude. If the 
maximum displacement amplitude is exceeded significant 
damage to the surround occurs. Once the desired amplum 
is achieved the refrigerator runs unattended until Sie 
State is reached. In this case, steady state is defined as 
less than a one to tuo percent deviation in the cold end 
temperature over a 90 minute period. This is determine EEEn 
evaluating the continuous computer monitor disploM 
temperatures, Which shows approximately 390 minutes of data 
history when set ap for display at four minute intel mins 
This is the basic procedure for startup and setup fo EE 
(applied) heat load data point. In general equilibration 
takes from 8 to 24 hours depending on the pressure 
amplitude setting. 

We have used two different methods for obtaining 
data with heat loads applied to the cold end. Them 
method is that used by Hofler (1986), where àa desi ns 
heat load is applied to the heater, at the cold end, and 
equili-bration is based on temperature stabilization, as 
discussed above for the no load case. This heat load is set 
by applying a fixed DC current and voltage to the ia 
coil. It was determined after tuo data runs that a Detter 
method might be available. This better method, designed and 


set up by Dr. Hofler, invloves a temperature feedback 
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ኺብ፡፡ያ፡ኸ፡፡ Network which allows setting a desired Steady state 
cold end temperature. The controller then establishes the 
current flou to the heater that is necessary to achieve and 
maintain that temperature. The equilibration point is no 
longer based on temperature but is instead based on the 
heater current deviation being less than approximately 
7@.5mA over a 9@ minute period. This corresponds to being 
Within better than ¿1% over the 90 minute period, for the 


currents we have measured. The control netuork involves a 


rather complex Fortran computer program, Written by Dr. 
Hofler, that reads the heater thermocouple, at an operator 
selected time interval, and then performs a crude time 
integral of the difference between the monitored 


memperature and the selected set point temperature. The 
program’s integrator output variable is used to numerically 
control the heater current from the HP6633A system D.C. 


power supply. The program allows the operator to select the 


control time interval (usually S seconds), the periodicity 
Of display (usually 4 minutes), the desired temperature, a 
proportional constant in volts/deg (usually 0.0), and an 
Mes Tataon Constant in volts/deg-sec (usually 0.01). The 


latter two of these control the stabilization time and the 
degree of stabilization. The usual values that were given 
۳۷۲۳ 1055-5 That Have Shown the best results during our 


experiments. 
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TRIS temperature feedback control method was 
determined to be better for the STAR project for two basic 
reasons. First, it provided a more rapid equilibration 
because the refrigerator was loaded more heavily at first, 
to drive the temperature of the cold end to the desired 
temperature, and then the load was let off only as 
necessary to maintain that temperature. This reduced the 
equilibration time to 6 to 8 hours in most cases. Secondly, 
this method conserves power due to the reduced equil- 
ibration time. This is considered to be vital to the STAR 
project due to its reliance on batteries for power |e 
gives it a limited operating lifetime. 

Jo Da tamkec ondas 

Once steady state has been established a data 56 iS 
recorded in a notebook. ል data set consists of nine 
readings, as well as date and time. First, the temperatures 
of the hot and cold heat exchangers are recorded. There are 
two thermocouples for each heat exchanger, one being Wau. 
center and one at the wall. Due to an intermittent failure 


^ 


01 the thermocouple at the center of the cold اا‎ 
exchanger we have only recorded the other three 
temperatures. Then we record the drive frequency, TEN 
pressure, the amplitude of the dynamic pressure signal and 
tHe amplitude of the lock-in velocity signal. Then, 


depending on which method of operation was used, we record 


the DC heater voltage and current, or just the DC heater 


56 


current. For the latter case we use the known heater 
resistance versus temperature characteristic to calculate 
the heater load (instead of voltage times current). After 
the data are recorded we change either the dynamic pressure 
amplitude or the heat load, and wait for the system to 
equilibrate. A Standard data run for a given gas type and 
stack involves 8 to 10 data points and takes about a week, 
from start to finish. The difference in the number of data 
points is based on the different temperatures achievable at 
different pressure amplitudes, While staying within the 
152013 * 5 ET ፡ amplitude limit. 

In addition to the above data, we also record the 
vacuum in the system and the temperature and resonant 
frequency of a gas analysis tube, similar to that described 
ከጠ ۵ و۱52۵‎ et. al. (1986), and Garrett, et. al. (1986), 
nue attached to the system. These are monitored just 
for completeness and are not directly involved in any 
Peta2ormance evaluation. With the exception of one data run, 
EN | our gases had analysis Certifications for purity, 
and these values uere used for calculations. In one case we 
mixed our own helium-argon gas and used the gas analyzer, 
1111م‎ by Ds Hofler, to determine E mixture 
concentration. 

The final bit of data that is recorded involves the 
۱۳۲ n Ch the engine warms up, after being fully cooled 


down and turned off. This data is analyzed and together 
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With the total heat capacity of the cold end of the SyS 40 
is used to determine a value for the heat leak out of the 
refrigerator. The heat leak determined by this method is 
then reduced by 2.4 mw/*c, which is the portion that has 
been determined to be due to the thermal conductance of the 


stack and gas, Which is internal and intrinsic to the 


This heat leak is then included in the total heat‏ ترما ره وميه 


load. 
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V. DATA ANALYSIS AND RESULTS 


In this chapter we will discuss the methods and equa- 
tions that were used in analyzing our data. This will be 
followed by a numerical and graphical presentation of 
various aspects of our experiment. It should be kept in 
mind that the purpose of these experiments was to improve 
the coefficient of performance (COP) of the prototype 


refrigerator, built by Hofler (1986). 


A. DATA ANALYSIS 

All of our data analysis was done through the use of a 
Microsoft Multiplan spreadsheet program and graphical 
display of the numbers that were thereby generated. This 
method was used because we had voluminous amounts of data 
where many of the required calculations were repetitive in 
nature. Taking into account all of the quantities, both 
recorded and calculated, we ended up with a spreadsheet of 
25 columns by 122 rows. It should be noted that many of 
these quantities were calculated for informational purposes 
and are not directly related to our final evaluation of the 
refrigerator’s coefficient of performance. 

Our analysis consists of comparing a normalized 


(normalized to the Carnot COP for the same temperature 
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difference) coefficient of performance, for the refrig- 
erator With modifications (COPR), to the same quantity for 
the original Hofler refrigerator. Since we have modified 
both the thermoacoustic stack and the gas type, we have 
three basic areas for comparison. We will compare the 
results with different stacks and the same gas type, with 
different gases and the same stack, and then with different 
stacks and gases combined. In order to make these 
comparisons we have created plots of COPR versus total heat 
load (Qtotal). In addition to these we have created plots 
of temperature ratio (Tc/Th) versus Qtotal, to validate our 
data. Based on the experiments of Hofler (1986), these 
latter plots are expected to be linear. 

In order to create these plots we need to know the 
values for COPR, Qtotal, Tc and Th. The values for Tc and 
Th are simply rau data taken during the experiments. We 
have used the temperature at the center of the hot heat 
exchanger as Th, and the temperature at the wall of the 
cold heat exchanger as Tc (due to the intermittent failure 
of the cold heat exchanger’s center thermocouple). 

The value for Qtotal has two components. The first is 
the column labelled QHEAT in Appendix C. This is the power 
applied to the heater, located below the cold heat 
exchanger, in watts. It is calculated by either multiplying 
the heater DC voltage times the heater DC current, or by 


multiplying the heater resistance by the square of the 
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heater DC current  (PsVI or P-I*R for DCO. The voltage and 
current values are raw data. The resistance values were 
taken from a plot of Tc vs. R created from the two sets of 
data where we recorded both heater voltage and current 
EVI. This plot is shown as Figure V-1. Both of these 
methods were used as discussed in Chapter IV. The second 
component of Qtotal is the heat leak, which is labelled 
HEATLK in Appendix C. As discussed in Chapter IV, this 
value is determined by analysis of the rate at which the 
refrigerator warms up once the power is turned off. Based 
on this warmup data and using the heat capacity of the cold 
portion of the resonator a value for the total heat leak to 
the cold end of the refrigerator is determined (in mW/*c). 
This number is then reduced by 2.4 mW/?c, uhich is the 
thermal conduction of the stack and the gas. The resulting 
quantity is the conduction external to the stack and is 
accuarate to Within 8%. The value for Qtotal is then the 
sum of QHEAT + HEATLK. The error in QHEAT is negligible, 
but the error in HEATLK is not. The error in Qtotal is 
about 8% when .QHEAT is zero and is about 1.95% to 2% for 
values of Qtotal near heat loads of optimum efficiency. 

The remaining quantity necessary for the calculation of 
COP is the acoustic power (or work). We use a relation 
Similar to that derived by Hofler (1986). The difference 
being that we no longer use the velocity signal from the 


bandpass filter; using the velocity signal from the lock-in 
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Figure V-1. Temperature vs. Resistance 
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instead. Also, we have added tuo correction terms to the 
equation. We are using an equation of the form: 


(1/Mp ° VpV, [C17 5,) 


(HH y) yy 


where Mp is the dynamic pressure transducer sensitivity, 


is the volume velocity sensitivity, Vp is the measured 


is the measured velocity signal, (i- 


W= 


HT 


pressure signal, VIT 


Ep? is an accelerometer non-linearity correction term, and 
SUV is a volume velocity calibration correction term. 
‹,=(7,4፳10 | Sec cmp a, where a is the acceleration in 
cm/sec? (x10 7), as shown in column 13 of Appendix C. 
Suy 79. 95 (Cf 0501 where f is the frequency in hertz, as 
shown in column 4 of Appendix C. The ratio 4 Hp is the 
result of the small cavity volume velocity calibration, 
discussed previously, and is therefore just a number that 
we put into the equation prior to calculation. Since the 
refrigerator was disassembled on several occasions during 
our experiments this calibration was performed several 
times. The calibration value changes over time as the 
driver is used and therefore this number Ay Hp is Corrected 
in our equations for all data taken after each calibration. 
The quantity Mp is also a calculated number that is input 
into the equations, since it is just the dynamic pressure 
transducer sensitivity. The actual equations that we have 
used are of the following form: 
5 


W-59810(1-7.4x10 a)V4(0.003V, ,)/(1839. 8-0. 05(£ -550)), 
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Where 598190 is C1/H,)? including unit conversions, (1- 
7.4x10 7a) is a correction term for accelerometer non- 
liearities, as discussed above, (0. 003V, 1) is the measured 
velocity signal (uith unit conversions in the coefficient), 
18898 15 ) ۲/1 ( the result of our volume velocity 
calibration, and (@.905(f-550)) is a correction to the 
volume velocity term due to the calibration being frequency 
dependent, as discussed above. The accelerometer correction 
(17 54) is about 2% at most, and the frequency correction 


(yy? is about 1% at most. The frequency correction was 
obtained for the case of pure helium. However, given the 
much larger frequency range of the gas mixtures, the amount 
of correction in some cases is too small, giving a 
calibration error of about 1%. The total worst case error 
for W is about 5% with 3.5% being a more typical value. 
Using the Qtotal and W values calculated by the above 
methods we get COP=Qtotal/W. To get COPR we need to divide 
this by the Carnot COP for the same temperature difference. 
The Carnot COP is known to be COP. arnot 1 0/ (Ih- Toe), Which 
means COPR can be derived as COPE=COP ((Th/Tc)-1). This is 
the equation we have used to determine COPR. On our 
Spreadsheet, in Appendix C this is the COPR2 column. The 
worst case error in the absolute values for COPR may be as 
high as 15% however, typical values are about 5%. Also, the 


relative error involved in comparing data taken under 


different conditions is probably considerably less than SX. 
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We now have all of the values necessary to evaluate the 
effects of our modifications on the refrigerator’s 


coefficient of performance. 


B. RESULTS OF EXPERIMENT 

ås stated earlier, we will break this discussion into 
three major sections. First we will discuss the results of 
our stack modifications. This will be followed by a 
discussion of performance with different gases and the same 
stack. We will conclude with the overall improvements using 
our best ርም combination. Also, a brief discussion 
concerning the results with different pressure amplitudes 
and different resonator lengths will be provided. 

Since we have data from two new stacks and two gas 
mixtures, as well as two different pressure amplitudes and 
two resonator lengths, available for comparisons, we must 
be selective in our presentation, So aS not to have an 
excessive number of plots that show the same result. To 
fulfill this goal we will present one or two sample sets of 
curves, with analysis, for each major category of interest, 
and provide numerical results for the other cases uithin 
that category, uhere apprpriate. In general a set of curves 
includes a plot of COPR versus Qtotal and a plot of Tc/Th 
versus Qtotal. 

Prior to presenting any of our figures ue shall intro- 


duce the notation that is used in the titles and legends. 
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All of our plots are designed to provide data comparisons 
and as such each contains more than one curve. The legend 
lists the specific curves using the following notation: 

GAS TYPE--STACK--RESONATOR LENGTH--PRESSURE AMPLITUDE, 
where abbreviations have been used in all cases. For gas 
types we have helium, helium-argon, and helium-xenon which 
are represented by He, HeáÁr, and HeXe, respectively. For 
the stacks ue have old, neu, and dual uhich are represented 
by O, N, and D, respectively. The resonator length is 
designated by s or 131 for short or long and the pressure 
amplitude is represented by the percent, relative to the 
mean pressure (Po/Pm), by the number 2 or 3. 

1. Stack modification results 
To illustrate the effects of our stack modifications 
we have elected to compare the data for pure helium, at 2% 
pressure amplitude, with the short resonator installed. 
Figure V-2 shows the temperature ratio (Tc/Th) versus total 


heat load (Qtotal) curves for all three stacks under these 


conditions. It is fairly clear from the curves that the 
data are essentially linear for all three cases, as 
expected. The curves also shou that for a given temperature 


ratio the old stack can support more heat load than the 
dual stack and the neu stack can support more heat load 
than either of the others. This statement can be reversed, 
such that for a given heat load the temperature spanned is 


greatest for the new stack and lowest of the dual stack. 
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Figure V-2. Temperature Ratio vs. Qtotal, For Helium, Short Resonator 
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These curves also show that our data are following the 
expected trends for the refrigerator, in general. Figure V- 
3 shows COPR versus Qtotal for all three stacks. It iS 
clear from this Figure that the COP is significantly 
improved with the new stack and vastly degraded with the 
dual stack. Evaluating the difference at the peaks of the 
curves, we see that the dual stack has a peak COPR of 0.034 
While the old stack has a peak at 0.067 and the new stack 
peaks at 0.078. This shows a 16% improvement for the new 
Stack and a 51% reduction for the dual stack. Conducting 
this same analysis for the remainder of our data sets 
yields an average improvement of 14% with the new stack and 
an average reduction of 57% with the dual stack. 

It is clear from these results that the dual Staki 
has not performed as predicted while the new stack has 
provided significant improvements. It is believed that the 
basic theory used to design the dual stack is valid, in-so- 
far as the concept is concerned. Houever, it is obvious 
that we have underdesigned the stack by either omitting an 
important factor or over simplifying. Át the time of design 
and construction of the stack, it was our belief that plate 


spacing, so long as it was not too small, would not have 


any detrimental effects. Based on this belief, and an 
initial evaluation of the new stack’s performance, we 


decided to build the dual stack with the hot end having the 


same fishing line size as that used in the new stack. 
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Figure V-3. COPR vs. Qtotal for Helium, Short Resonator 
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However, with the dual layer at the hot end and the thick 
fishing line at the cold end, we ended up with a 20% 
increase in plate spacing at the hot end and a 163% 
increase in plate spacing at the cold end (relative to the 
spacing in the new stack). Based on the extremely poor 
results with the dual stack it appears that the magnitude 
of the plate spacing is more critical than we had 
originally thought. The theory, presented in Chapter II, 
shows that a plate spacing of 2 to 4 thermal penetration 
depths, with the entire resonator cross-section filled with 
plates 15 desired for the best heat pumping 
characteristics. However, our analysis of the theory did 
not predict this massive drop in COP with a larger plate 
spacing. To some extent this plate spacing effect is 
illustrated by the new stack’s performance, as well. 
Although our intuitive reasoning for the design of the new 
stack predicted improved performance, the results are 
better than expected. Combining this with the extremely 
poor performance of the dual stack, and the fact that the 
new stack has the smallest plate spacing of the three 
Stacks, we can see that plate spacing iS much more 
important than previously thought. Rough calculations for 
helium show that the average plate spacing, in number of 
thermal penetration depths is 3.3 for the new stack, 3.5 
for the old stack, 4.0 for the hot end of the dual stack, 


and 8.8 for the cold end of the dual stack. 
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2. Results for Gas Mixtures 

For this analysis we have elected to present two 
sets of data, both of which are for the new stack. This is 
primarily due to the fact that the new stack is the only 
one for which we have any helium-argon data. Also, based on 
the stack analysis, the new stack is the best one 
available. For this analysis we will discuss the results 
for 2% pressure amplitude and both the short and long 
resonator. The long resonator consists of a one inch spacer 
that lengthens the resonator and thereby changes the 
ሽክ On of - the stack realtive to the acoustic standing 
wave. 

Figures V-4 and V-95 Show the plots of Te/Th versus 
EN ለጠ TOL the short and long resonators, respectively. It 


is clear from these curves that all of the data are linear, 


as expected. It is interesting to note that in Figure V-2 
all of the curyes for pure helium were essentially 
parallel, even though the data was for different stacks. 
Now, in Figures  V—-4 and V-5 ue see that the curves are no 


longer parallel for the different gases, even though the 
data are all for the same stack. This is believed to be 
related to the fact that the heat pumping capacity is 
related to the PU product (pressure times volume velocity). 
For helium gas at constant pressure amplitude the PU 
product will vary in a consistent manner, with respect to 


temperature, regardless of stack design. For different 


||| 


1۱/7 
MIA 
HUOT AI 
wl LTE 
1 
LLLA LLL LLLI gs: 
|||/ሪ/|||||| |  ሙ 
“ዘ ተገገገገገተ 


Y 





0.70 


7 


ALL i ttt tt 
ot LLL TLL] 


Q total (watts) 


Figure V-4. Temperature Ratio vs. Qtotal, New Stack, Short Resonator, All Gases 
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Figure V-5. Temprtature Ratio vs. Qtotal, New Stack, Long Resonator, All Gases 
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gases however, the velocities will vary because of 
different sound speeds, and therefore even at constant 
pressure amplitude the PU product will vary differently, 
producing curves with different slopes. 

The plots of COPR versus Qtotal for the short and 
lone resonator are shown in Figures V-6 and ۷-7 
respectively. It is clear from all of these curves that 
helium-argon produces a certain degree of improvement in 
COP, over pure helium, and that helium-xenon produces the 
highest COPR. For the short resonator curves of Figure V-6 
we see that the peak COPR values are approximately 0.078 
for helium, 0.102 for helium-argon, and 0.129 for helium- 
xenon. From these values we see that helium-argon gives a 
31% increase over helium, while helium-xenon provides a 65% 
improvement over helium and a 27% improvement over helium- 
argon. It should be noted that the peak COPR values are 
shifted to lower heat load values for the gas mixtures. 
This is believed to be due to a reduction in this PU 
product that was discussed above. Similar calculations for 
Other stacks and other pressure amplitudes show an average 
improvement in the peak COPR values of 61%, when comparing 
helium-xenon to helium. There is no other short resonator 
data available for helium-argon. The long resonator curves 
of Figure V-7 have been included primarily because they 
illustrate the peak COPR value that has been measured to 


date of 90.197 for helium-xenon. Although the helium-argon 


|| + 


a ۲۳۵ ۸! و‎ 2 

o HexeNs2 

a ۲6۸۵۲ ۸! و‎ 2 
0.0 0.5 1.0 1.5 2.0 25 30 


Q total (watts) 


e 
nO 


o 


0 


ን 


0.0 





O) 


0.0 






Figure V-6. COPR vs. Qtotal,New Stack, Short Resonator, All Gases 
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Figure V-7. COPR vs. Qtotal, New Stack, Long Resonator, All Gases 
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curve appears someuhat incomplete ue feel safe in assuming, 
based on all of our other curves and data, that the curve 
is within a few percentage points of its peak value at the 
last plotted data point. We will therefore use that datum 
for our peak value. These curves show the same general 
trends as those of Figure V-6, however the magnitudes of 
the improvements are somewhat reduced. For peak values of 
0.132 for helium, 9.1595 (from data of Appendix C) for 
helium-argon and 9.197 for helium-xenon, ue see that helium- 
argon only produces a 21% improvement over helium and 
helium-xenon produces a 49% improvement over helium with a 
24% improvement over helium-argon. It is believed that the 
improvement percentages With the long rresonator are 
reduced because the long resonator itself produces about a 
62% improvement over the short resonator, thereby placing 
the refrigerator closer to a maximum possible or limiting 
COPR value. Thus the absolute improvements ue achieved have 
00 وا تا دا‎ ۳161 ۷۵086 wise. 
3. Overall Improvement 

In this section we present the overall result of the 
thesis. That is, the overall improvement in the COP of the 
prototype refrigerator, taking into account both the stack 
modifications and the binary gas mixtures. To illustrate 
this case we have elected to compare the results with the 
Old stack and pure helium at 2% pressure amplitude to those 


achieved with the new stack and helium-xenon at 2% pressure 


|| 7 


amplitude. These sets of data were selected because the 
original refrigerator used only the old stack and pure 
helium, while our best results were achieved with the new 
stack and helium-xenon. We have not used the long resonator 
data, even though that gives our highest COPR, because we 
did not take any long resonator data with the old stack. 
However, for a number for the old stack, long resonator, 
pure helium setup we can use the 9.126 peak COPR given by 
Hofler(1986) in his doctoral dissertation. Comparing this 
to our best case peak COPR of 0.197 we get a 56% 
improvement. 

For the comparison of our experimental data we will 
present the COPR versus Qtotal plot only, since the Tc/Th 
versus Qtotal plots have already been discussed for this 
data. Figure V-8 shows our plot for overall improvement. 
With peak COPR values of 0.129 for the helium-xenon and 
0.067 for the helium ue see that the overall improvement in 
the refrigerator’s coefficient of performance, due to our 
modifications, is 93%. 

4. Other Effects of Interest 

There are two other areas of interest that have been 
mentioned briefly during our discussions of results. These 
are the effects produced by changing resonator length and 
those produced by changing the pressure amplitude. 

In Figure V-9 we present the data for all three 


gases in the new stack with both short and long resonators. 
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Figure V-8. COPR vs. Qtotal, Showing Overall Performance Improvement 
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Figure V-9. COPR vs. Qtotal, Showing Resonator Effects 
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From this We can see that just the shift from short to long 
resonator produces increases in GCOPR as follous: 70% for 
helium, 59% for helium-argon and 52% for helium-xenon. This 
rather large increase in COP is due primarily to the shift 
in the position of the stack relative to the Standing wave. 
With the long resonator installed, the hot end of the stack 
is moved closer to the velocity antinode and the heat 
transfer power is increased Significantly, thereby 
producing the large increase in the COP and the observed 
shift of the COPR peak towards higher heat loads. 

In Figures V-10 and V-11 we present a comparison 
between both gas types and pressure amplitudes. In Figure V- 
10 we have helium-xenon with the short resonator, both 
pressure amplitudes and all three stacks. In Figure V-11 we 
have all of the same curves for pure helium. In all of 
these curves it is obvious that at higher pressure 
amplitudes the peak COPR value has shifted toward higher 
heat loads. Based on this and the discussion of Figures V-4 
and V-5 it is clearly shown that heat transfer power (or 
heat pumping capacity) increases with pressure amplitude 
(due to the PU product increasing). It is interesting to 
note that the peak COPR value is in fact lower, at T 
pressure amplitudes, for all of the cases for uhich ue have 
data, With the exception of the pure helium, short 
resonator ad stack case. In his doctoral dissertation, 


Hofler reported that the peak COPR value increased with 
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Figure V-10. COPR vs. Qtotal for HeXe, Short Resonator, All Stacks and Amplitudes 
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Figure V-11. COPR vs. Qtotal for He, Short Resonator, All Stacks and Amplitudes 
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pressure amplitude for all of his data. In fact our damn 
With his stack and helium also follow this trend, 
Additionally, Hofler (1986) reported: "there is also some 
indication that the magnitude of these maxima approach some 
constant value at higher amplitudes." Our results for gas 
mixtures and other stack designs tend to contradict this, 
and lead toward a rather alarming conclusion that there is 
some relatively low pressure amplitude at which the COPR 
maxima peak and above which the COPR maxima drop off again. 
This Seems to indicate that there may be "Some 8 
linearity, such as turbulence, that limits the best COP's 


to relatively low amplitudes. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The primary purpose of this thesis iS to provide 
modifications for the Space Thermoacoustic Refrigerator 
(STAR) that will improve its thermodynamic efficiency, 
specifically its coefficient of performance (COP). Through 
the use of a helium-xenon binary gas mixture and a new 
stack design, with less string and reduced plate spacing, 
we have produced a 93% increase in the peak value of COP. 

Along with this improvement we have raised Several 
questions that are as yet unanswered. Based on our results 
with a long resonator, relative to those with a short 
resonator, it would appear that there iS Some optimum 
resonator length and design that must be found. Further 
testing with different resonator Spacer thicknesses, must 
be performed to resolve this issue. The unexplained drop 
in peak COP relative to Carnot (COPR), with increasing 
pressure amplitude opens the door for more testing with 
different setups and many pressure amplitudes. It is clear 
that the binary gas mixture of helium-xenon is superior 
under all of the conditions that we have observed. However, 
the failure of the dual stack and the better than expected 
results with the new stack leave the final stack design as 
an open item. It is clear that another dual stack and a 


different version of the new stack should both be built, 
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with smaller plate spacings, to determine the optimum stack 
design. It is the answers to these question that will yield 
the design capable of producing the maximum thermodynamic 
efficiency for the Space Thermoacoustic Refrigerator 
(STAR). 

Although this project deals specifically with the space 
application of our refrigerator, there iS no reason that it 
should be limited to that application. Due to the lack of 
moving parts that are present in most other cooling engines 
(i.e. sliding seals) this unit should be extremely 
reliable. It also has reasonable efficiency relative to 
most commercial refrigeration devices. The fact that it 
does not use 581657 1 አዘዘች ለይ is a great advantage, 
especially With the close scrutiny that this ozone 
destroying substance is currently receiving. Therefore, 
this type of refrigeration device could, and should be 
developed for commercial use and should be considered as a 


possible candidate for virtually any type of cooling need. 
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APPENDIX A.  PRANDTL NUMBER PROGRAM 


This appendix provides a program and documentation for 
the calculation of Prandtl numbers for binary gas mixtures. 
The progran, Written by the author, is designed to 
calculate Prandtl numbers for all concentrations of a 
binary mixture of Helium and one other gas. The program 
is written in Fortran for use on the Naval Postgraduate 
school mainframe computer, with a Watfor-77 compiler. Two 
Output files are created when the program is run. The 
first, called "PRANTL OUTPUT Al", is a listing of Prandtl 
numbers Versus helium" fraction for further use in a 
plotting program, to create Prandtl number versus helium 
fraction curves for the gas mixture selected. The second, 
called “PRANTL RESULTS Al" is a tabular listing of Prandtl 
number, thermal conductivity (in cal/cm-sec-K), and dynamic 
viscosity (in gm/cm-sec) versus helium fraction for the 
selected gas mixture. The listings go from pure helium to 
zero percent helium in one percent intervals to provide 
easy access to any specific mixture. All of the theory and 
equations necessary to create the program, other than 
knoWledge of a programming language, can be found in 
"MOLECULAR THEORY of GASES and LIQUIDS" by Hirschfelder, 
Curtiss, and Bird, John Miley & Sons, Inc., New York, 1954. 


The primary source was chapter eight of the text, hereafter 


7 


referred to as Hirschfelder, although section and equation 
numbers Will be Specified where necessary for 
Clarification. 

ÀS Written, this program uas designed for personal use, 
exclusively to derive the data necessary to create the 
previously mentioned curves of Prandtl number versus helium 
fraction. It has been modified to include the second 
tabular output for completeness and further use by others 
in this field. The program is currently written to require 
hard input and therefore requires editing of the main 
program for each change of gas type or temperature. In 
addition to satisfying the specified needs, this method was 
used due to there being thirteen input parameters that vary 
With gas type and temperature, five of which require table 
lookups in Hirschfelder. This program could be modified for 
incorporation into a larger program as a subroutine, 
provided the specific gas types and temperatures of 
interest are completely specified and all of the input 
parameters are looked up and incorporated into a data array 
for reading. This has not been done to date as the current 
useage did not require these features. 

As previously stated the primary source for this progran 
is Chapter eight of Hirschfelder, which provides the 
equations necessary +0 calculate Various transport 
coefficients, such as the coefficient of viscosity and the 


coefficient of thermal conductivity. The equations 
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presented in chapter eight rely heavily on the elementary 
theory of transport phenomena presented in section 1.2 of 
the text as well as the rigorous kinetic theory for 
monatomic gases presented in chapter seven. The intricate 
details of this theory will not be presented here, however 
a discussion of the methods used, including the main 
assumptions and restrictions, will be provided as well as a 
definition of each of the terms or symbols in the program. 
The development of the kinetic theory for gases is based 
on knowledge of the distribution function that represents 
the number of molecules of a specific species, in a unit 
volume element about a point in space, with velocities ina 
unit range about a specified velocity, at a given instant 
in time. Since the usual definitions for transport 
coefficients only apply under conditions of equilibrium (or 
only slightly different from equilibrium) this restriction 
applies to the results. Under this limit the distribution 
function 15 nearly Maxwellian and is solvable by a 
perturbation method developed by Chapman and Enskog. The 
solutions are then used to obtain expressions for transport 
coefficients in terms of a set of integrals, Omega(l,s), 
which involve explicitly the dynamics of molecular 
encounters. The Chapman-Enskog theory relies on several 
assumptions and therefore has limited applicability. Each 
assumption will be briefly discussed and related to the 


specific situation of our experiment. 
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Only binary collisions are considered, Which limits us 
to densities lou enough that three body collisions are 
unimportant. The theory predicts that viscosity is 
independent of pressure at constant temperature. 
Experiments on Nitrogen have shown that only a 4% deviation 
in viscosity occurs as pressure is varied form 1 to 60 
atmospheres. Our operating condition of 1@ atmospheres is 
well within this and should therefore closely follow the 
theory in this respect. 

The theory further limits itself to Classical Mechanics 
and thereby precludes low temperature situations where 
quantum effects become Significant. These quantum effects 
are less than 1% for helium above 200 Kelvin (and even 
smaller for heavier isotopes). The coldest temperature we 
have achieved to date has been 193 Kelvin under a no load 
condition so our normal operations are all above 200 Kelvin 
and we are not concerned with quantum effects. 

The Chapman-Enskog theory uses a series approximation, 
to the above mentioned distribution function, in solving 


the Boltzmann equation and then uses only the first 


approximation, which is only valid When the gradients of 
the physical quantities (i.e. density, velocity, 
temperature) are small. The changes of these properties 


over one mean free path must be small with respect to 
unity. Based on the mean free paths for gases at pressures 


above one atmosphere this condition is satisfied except 
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under conditions of extreme gradients such as those present 
in shock waves. Since we do not have any shock waves we are 
not restricted by this. 

The theory further assumes that the dimensions of the 
containing vessel and any obstacles within it are large 
compared to the mean free path, which provides no 
limitations for our gases and our pressures. Although the 
Chapman-Enskog theory applies strictly to monatomic gases 
the results have been shown to be good, even for polyatomic 
gases, provided that the molecules are not to non- 
spherical. Finally, the development of the theory results 
in the set of integrals Omega(l,s), which are present in 
all of the equations for transport coefficients. These 
integrals, in turn are dependent upon the potential 
function that represents the molecular interactions. 
Throughout the derivation in Hirschfelder, and for all of 
the — o used here, the Lennard-Jones (6-12) potential 
is assumed. 

Since some of the other terms and symbols depend upon 
the Omega terms in the program, we will start with a 
discussion of where the Omega terms come from. These 
Omega(l,s) terms represent a set of integrals, linear 
combinations of which are used to evaluate the more complex 
bracket integrals arrived at in solving for the expansion 
coefficients of the Sonine polynomials. The transport 


coefficients are expressed in terms of the Sonine 


|.3| 


polynomial expansion coefficients and can then be expressed 
in terms of these Omega(l,s) integrals. In Hirschfelder 
there is also an Omega-star(l,s) symbol Which represents 
the deviation of a particualr molecular model (in our case 
the Lennard-Jones (6-12) potential) from the idealized 
rigid sphere model of molecular interactions. All of the 
Omegas in the program represent an Omega-star value. The 
Omega integrals are functions of the reduced temperature 
(Tx=kT/¢) and are tabulated in table I-M of Hirschfelder. 
Where more than one variable of a given type exists, 
such as Sigmal, sigma2, and sigmi2, the 1 and 2 represent 
the value of the quantity for pure gas 1 or pure gas 2 
(which is always helium for our purposes), respectively. 
The 12 represents the value for gases 1 and 2 together. 
This is true for all cases except ETA12 and LAM12. These 
two variables represent artificial quantities for a 
hypothetical pure substance with molecules of molecular 
weight 2(M1)(M2)/(M1*4M2) that interact according to a 
potential curve specified by interaction parameters sigmal2 
and የ12. These parameters are defined in terms of the 
Lennard-Jones (5 — 120) potential in Chapter one of 
Hirschfelder, Where sigma is the distance of A 
approach (in angstroms) of tuo molecules uhich collide with 
zero initial relative kinetic energy and * is the maximum 
energy of attraction of the two molecules. Table I-A of 


Hirschfelder provides values for Sigma and ¢/k for various 
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substances, where k is the Boltzmann constant and ¢/k has 
units of Kelvin. The parameters f12 and sigm12 are arrived 
at through empirical combining laws which relate force 
constants between unlike molecules to those between like 
molecules. These laws are sigm12=(1/2)(sigmal+sigma2) and 
€12=( +425 The quantities A12 and Bi2 are merely 
frequently  encounterd ratios of different Omegas that have 
been tabulated separately as functions of Tx in table I-N 
of Hirschfelder. 

The remaining terms or symbols in the program are fairly 
Straightforward, based on the above explanations and each 
type will now be defined: 

PR=prandtl number 
HE=helium fraction 
LAM-thermal conductivity coefficient (in cal/cm-sec-K) 
ETA=coefficient of dynamic viscosity (in gm/cm-sec) 
M=molecular weight 
CP=specific heat capacity at constant pressure 
(in cal/gm-K) 
T=absolute temperature (in Kelvin) 
Al=fraction of non-helium gas 
AZ=helium fraction 
T=counting variable for the Do loop 

The output parameters are Prandtl number, helium 

fraction, thermal conductivity, and dynamic viscosity. They 


are represented by PR(I),HE(I),LAMIX(I), and ETAMIX(I) 


53 


respectively. The remaining parameters are all merely 
mathematically defined intermediate steps in the 
calculations of the output parameters and have no physical 
meaning. The terms used for viscosity are XETA, YETA, and 
ZETA, and the terms used for thermal conductivity are Ul, 
U2, UY, UZ, XLAM, YLAM,; and ZLAM. 

This section will now provide the most rececnt version 
of the program. Since the explanation of the input is 
fairly lengthy it will be covered first, rather than 
disrupt the flow of the program. This will include a list 
of all required input parameters, their sources where 
appropriate, and a sample set of inputs that were used. 
There are thirteen input quantities, A1l2, B12, OMEG12, M1, 
M2, CP1, CP2, SIGMA1, SIGMAZ, SIGH1Z, OMEGA], OMEGAZ, AND 
T. After selecting Mi, M2, and T you have specified your 
gas combination and temperature of interest. Values for CP 
may be obtained from various sources, such as the CRC 
tables, thermodynamics textbooks, or standard tables of 
thermodynamic properties. The values for SIGMA as well as 
¿/k are obtained from table [I-A of Hirschfelder. Using */k 
for each gas and £12=( €1xé2)” you calculate T1x, .و12‎ and 
[1 12 (where 3 Tx-kT/:). Now i the appropriate Tx you go 
to table I-M of Hirschfelder for the OMEGA values and table 
I-N for Al2 and B12. This provides a complete set of input 
quantities. The following table is a sample set of inputs 


for two gas combinations at two different temperatures: 
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GAS COMBINATION AND TEMPERATURE 


He-Ar 


300 


39.948 


4.0026 


3.418 


6 


124 


10.22 


1.104 


0. 7485 


©. 12477 


1251-2 


35. 60 


8.427 


1.108 


1.0905 


0.8470 


the array variables are dimensioned to allou 


He-Ar 


200 


39.948 


4.0026 


3.418 


2.5936 


124 


10. 22 


TS 


0. 7467 


©. 12477 


1 2 


35 . 0 


5.8 


1. 2 


1.091 


0.3080 


He-Xe 


300 


1313 


4.0026 


4.055 


2556 


229 


10.22 


1.394 


0.7489 


0.0382 


12512 


48. 38 


6.201 


1. 1034 


1.0904 


0.8316 


He-Xe 


200 


13173 


4.0026 


4.055 


2:396 


229 


10.22 


1.700 


0.7467 


0.0382 


2312 


48.38 


4. 4 


1.098 


1.093 


0. 2 


LoD 


QUANTITY 


OF INTEREST 


M1 


M2 


SIGMA1 


SIGHAZ 


(£/K)1 


2 ( 1ع )2 


OMEGA 1 


OMEGA2 


CP1 


CP2 


IZ KO 


1 12 


A12 


B12 


OMEG 12 


PROGRAM PRANTL 


First 


C 


C 


C 


calculating and plotting the desired quantities. Then 
the real quantities are defined since some begin with 


letters that uould make them integers by default. 


DIMENSION ۳ ) 150۱ , ۲1۳ ) 150 ( , LAMIX( 150, ETAMIX( 150) 
REAL H1,HM2,0MEGA1, OHEGA2, OMEG12, LAM1, LAM2, LAM 12, PR, HE 


REAL LAMIX,ETAMIX 


Nou comes the input section where the counting variable 


is initialized to zero and all required data is input. 


A12=1. 1080 
8512-1. 5 


ONEG 12=0. 84790 


1-6 
M1=39.948 
M2=4.0026 


CP1=0. 12477 

CORZaden2 58.2 

SIGMA 1=3. 418 

6 .1182-2 1 5ت 

5 1 11 12-0 . 5 ) 5 ] 6۳۸ 1 +5 12( 


OMEGA 1=1. 104 
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OMEGA2=0. 7485 


T=300 


Now the output files are created and opened and headings 


are placed where appropriate. 


OPEN(67,FILE=*PRANTL OUTPUT A1?) 
OPEN(68,FILE-'PRANTL RESULTS A1?) 
WRITE(68,X)' THIS OUTPUT FILE IS FOR A MIXTURE OF 
XGASES WITH? 
۷ ۴ 1 15 ) 658 , 10111 , 2 
1 FORMAT(* ATOMIC WEIGHTS OF ',F10.7," AND”*,F10.7,” ,AND 
xk A °) 
WRITE(68,2)T 
CEUURMATCUDEBPERATURE OF ۰,۲ 11. ۲۰ DEGREES KELVIN.?) 
WRITE(68,*) 
WRITE(68, x)’ PRANTL HE THERMAL 
* DYNAMIC" 
WRITE(68, *) ” NUMBER FRACTION CONDUCTIVITY 
x VISCOSITY’ 
WRITE(68, x) 
WRITE(68, x)’ CAL/CM-SEC-K 
* GHM/CH-SEC' 


WRITE(68, x) 
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Nou the Do loop is opened to perform the required 
calculations. The equations are from chapter eight of 
Hirschfelder and the equation numbers are as follous: 
ETA1 and ETA2 equation 8.2-18 

ETA12 equation 8.2-21 

LAM1 and LAM2 equation 8.2-31 

LAM12 equation 8.2-35 

U1,U2,UY,UZ, XLAM, YLAM, and ZLAM equation 8.2-36 

XETA, YETA, and ZETA EQUATION 8. 2-2 


The equation for Prandtl number is ETAMIXXCPMIX/LAMIX. 


DO 10 X1=0.0,1.005,0.01 

1=1+1 
ETA1z2.6693E-S5XxSQRT(MIXTO)/(SIGMA1XSIGMA1XOMEGA1) 
ETA2z2.G6693E-S5X5QRT ( MZXT)/(SIGMAZXSIGMAZXOMEGA2) 
ETA 12=2. 669TE-OXSQRT (2xM1xM2x*xT/(M1+M2))/(SIGM12x 
*SIGM12xO0MEG12) 
LAM1=1.9891E-4xSQRT(T/M1)/(SIGMA1xXSIGMA 1x0OMEGA 1) 
LAM2z-21.9891E-AXSQRT(T/M2)/(SIGMAZXSIGMAZXOMEGA2) 
LAM12=1.39891E-4xSQRT (Tx(M1+M2)/(2x*xM1xM2))/((SIGM12 
XXX*2)x0MEG 12) 
Ul=(4xA12/15)-1(((B12/5)+(1/12))x(M1/M2))+(((M1-M2) 


XXX2)/ (ZXM1XM2)) 
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U2=(4xå12/15)-(((B12/5)+(1/12))x(M2/M1))+(((M2-M1) 
xxx2)/(2xM1xĦM2)) 

UYz((AXA12/ 150 XCCCM1- H2) XX2)/ CAXM 1XM2) 2) XCLAM 12ZXLAM 
۷ 127 ) ۸۵۸۲] 1۴1,۸۲2 ( ( ( - ) ) 12۸/5۱ + )1۸/12((- ۲ ) ) 12 127/32۷۸12 
#۷ ( - ) 257۸ ) 32*۸ 12( ( (7) ) )۳1-۳2 ( ۷2 (7 ) ۳۲ 1:۳2 ( ( ( 

UZ=((4xA412/159)x(((M1+H2)xx0)/(4xM1xM2))x((LAM12/ 
LANA) +(LAM12/LAM2))-1))-(B12/5)-(1/12) 

X2=1-X1 

XETA=(X1xXX1/ETA1)+(2xX1xX2/ETA12)+(X2xX2/ETA2) 

YETA=(3xA 12/9)x((X1xX1xM1/(M2xETA1))+(X2xxX2xH2/(Mi1x 
kETAZ))+(2xX1xX2x( (M1+M2)xx2)x(ETA12xETA12)/(4xM1xHM2x 
*ETA LXETAZXETA12))) 

ZETA=(3xA12/5)x( (X1xX1*xM1/M2)+(X2*xX2*xHM2/M1)+( (2x*xX 1x 
AXZIX(((((M1+M2)xx2)/(4xM1xM2))x((ETA12/ETA1)+(ETA12 
X/ETA2)))-1))) 

XLAM=(X1xXX1/LAM1)+(2xX1xX2/LAM12)+(X2xX2/LAM2) 

ሃፒልሽ= (፪ 1#፪ 1#ሀ 1/ፒልሽ 1) + (2#8 1#፳2#ሀኘ/ፒልቨ 12) + (፤2#8፪2#ሀ2/ 
XLAM2) 

ZLAM=(X1xX1*xU1)+(2xX 1*kX2x*U2Z )+ (X2*xX2xU2) 

PR(1I)=((14+ZETA)x(XLAM+YLAM)x( (X1*xM1xCP1)+(X2xM2xCPZ2) 
X))/((XETA+YETA)x(1+ZLAM)x((X1xM1)+(X2xM2))) 

HE(1)=X2 

LAMIX(1)=(1+Z2LAM)/(XLAM+YLAM) 


ETAMIX(1)=(14+ZETA)/(XETA+YETA) 


3 


Now we will write to the output files the data that has 


been calculated. 


WRITE(67,5)PR (1) ABI) 

WRITE (68,6)PR(1),HE(1), LAMIX(1),ETAMIX(I) 
5 FORMAT(F10.7,F10.2) 
6 FORMAT (ሇብ ዖን. 7, 12.2575 7 
19 CONTINUE 

STOP 


END 
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APPENDIX B. PLOTTING PROGRAM 


This appendix provides a program, written by the 
author, that will draw curves of Prandtl number versus 
helium fraction from the output file, “PRANTL OUTPUT Al", 
created by running the program "PRANDTL FORTRAN Ai" that 
was presented in Appendix A. This program must be compiled 
in FORTVS and is specifically designed to use DISSPLA on 
the Naval Postgraduate School mainframe computer. Various 
DISSPLA subroutines are called by the program and each 
will be explained With comment statements. DISSPLA 
requires the data that is to be plotted to be in an array 
format. Each of the arrays used here are set to a maximum 
of 150 points. If modifications are made to "PRANDTL 
FORTRAN Ai" such that more than 150 points are calculated 
the array sizes must be redefined. Two character strings 
are used in this program. They are NEW1$ and NEW2$, both 
of which are used to rename files. The first changes the 
input filename into the format required for DISSPLA and 
the second changes the file back to its original name for 
storage. The only parameters used are PR(I) and HE(I) 
which represent Prandtl number and helium fraction, 
respectively. Both of these quantities are real numbers 
and are defined as Such in the program. This completes the 


introduction and explains the first three lines of the 
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program. The remainder of this appendix uill be a listing 
of the program uith explanatory comments. Latest update of 


program, March 1988. 


Program GRFPRN 
characterx80 NEW1$,NEW2S 
dimension PR(150),HE(150) 


real PR,HE 


C The next three lines rename the input file to the 


C format necessary for DISSPLA and open it for use. 


NEW1$=’rename '//? PRANTL OUTPUT A1"//? file 
KXYZ1 A1? 
call excms(NEW1$) 


open(67,file=’XYZ1°) 


c In the next seven lines a Do loop is used to read 

c all of the data from the input file into the 

C Specified arrays. Then the input file is closed and 
C renamed it to its original name. The format line 

C ensures the data is read from the file in the same 


C format aS it Was written into the file. 


1 format(f10.7,f10.2) 


ኋጋ 152 1=1, 1525 
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read(67, 1)PR(I),HE(I) 
10 continue 
20 close(67) 
New2$=*rename °//’ file XYZ1 Al’//’ PRANTL 
xXOUTPUT A1’ 


call excms(NEW2$) 


The remainder of the program calls DISSPLA subroutines 
to create and plot the desired curves. Call tek618 

Sets up the required machine for the plotting and 

links you to the plotter. You must be on a tek618 
terminal to use this program. Call blouup is used to 
magnify the plot to near full page size, for clarity. 
Call page sets up the page siZe, in inches. Call area2d 
sets up a two dimensional plot area of the specified 
size, in inches. Call xname and yname allow you to 
label the axes. Call headin allows naming the plot. 
Call graf sets up the physical origin and the maxima 
for the axes. The ’scale’ portion of this call allows 
DISSPLA to select appropriate increments for the axes 
based on the minima and maxima specified. Call thkfrm 
and call frame place frames around the inner plot area 
and the entire plot, of a specified thickness. Call 
curve tells the plotter to plot the specified arrays of 
data points. The arrays must be called in the proper 


order (y-axis, x-axis) to get the correct curve. Call 


143 


c endpl ends the specific plot and sets up a new page for 
c another plot. Call donepl tells the device you are done 


C ات‎ ۲ 


Call tek618 

call blowup(1.5) 

call page(11,8.5) 

call area2d(8.,6.) 

call xname('helium fraction$',100) 

call yname(’Prandtl number$’, 1900) 

call headin(’Prandtl number vs. helium fractions’, 
*100,1.,1) 

call &raf(0.0, "scales. 190.0. d DESEE? 007) 
call thkfrm(.02) 

call frame 

call curve(HE,PR, 190,0) 

call endpl(9O9) 

call donepl 

Stop 


end 
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APPENDIX C. RAW DATA 


This Appendix provides a listing of our spreadsheet, 
discussed in Chapter V, including all recorded data and 
calculated values that we have used. Many of the values are 
not related to our coefficient of performance evaluation, 
but were calculated for informational purposes. 

The Appendix is broken down into three sets of three 
pages each, in order to fit the 24 column by 122 row spread- 
Sheet within the desired margins. Each set of three pages 
covers the entire 24 columns for the specified data sets, 
With nine columns per page. We will discuss each of the 
columns, omitting the first which specifies the data set. 

Columns 1,2,3,4,6,8,9,10 and 11 are all raw data and are 
self explanatory. Column 5 is the measured mean pressure of 
column 4, after conversion from milli-volts to psia (at 1.5 
psia per mv). Column 7 converts the measured dynamic 
pressure signal of column 6 to percent pressure amplitude 
(pg/ P,71586xC6G/C5). Column 12 is an acceleration term that 
is calculated for use in the uork equation as a correction 


for accelerometer non-linearities (acc=@.2848xC3xC8/170.5). 


Column 13 is the absolute temperature ratio 
(To/Ty=(C2+273. 15)/(C1+273. 15)). Column 14 is the applied 
heater load in watts (QHEAT=C9xC 10x0.001 or 


QHEAT-C10xC10xC117109, as appropriate). Column 15 is the 
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work in watts (W-59810x(1-7.4x10 730 12 (06:60 ۰ 003087] 1839. 8- 
0 . 05 )03-550( 1 as discussed in Chapter V. Column 16 is the 
coefficient of performance, neglecting the heat leak from 
the cold end of the refrigerator (COP=C14/C15). Column 17 
is the coefficient of performance relative to Carnot, 
neglecting the heat leak (COPR=C16((1/C013)-1)). Column 18 
is the heat leak in watts, as discussed in Chapters IV and 


V  (HEATLK=0.0064(C1-C2)). Column 19 is the coefficient of 


performance, taking into account the heat leak 
(COP2=(C 14+C 18)/C15). Column 20 is the coefficient of 
performance relative to Carnot, taking into account the 


heat leak (COPR2=C19((1/C13)-1)). Column 21 is the total 
heat load in watts (Qtot=C14+C18). Column 22 is the overall 


temperature difference (DELTA T=C1-C2). Column 23 is an 


approximate viscous penetration depth for the given 
temperature, pressure, frequency and gas type and an 
approximate value for dynamic viscosity (VISC 


PEN=SQRT ( (C1+273. 15) (0. 005014692852) /(C3xC5)), where the 
, coefficient takes into account all constants and the value 
used for dynamic viscosity). Column Z4 is an approximate 
thermal penetration depth calculated from & = 7% (THERM 
PEN-1.219820636xC23). The equations for viscous and thermal 
penetration depths use a single value for and ¢ within 


each data set, respectively, instead of a value calculated 


for each data point. This is done for ease of calculation. 


| +6 


The resulting values are therefore not precise but merely 
reasonable approximations. 

The spreadsheet of raw data and calculated values is 
provided on the next nine pages, and does not follow on 


this page for formatting reasons. 
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VHEAT 

V 

0,000 
0. 000 
6.870 
7.714 
11.897 
15. 780 
13.780 
15.780 
0. 000 
0.000 


0,000 
5.871 
9.706 
11.904 
8.416 
0. 000 
0. 000 
11.898 


VLOCKIN 
aVrms 


رمد 
241.6 


257.1 


277.0 


289.1 


307.0 
309.0 
3215.3 
2200 
526.5 


ATY 


ፈጋ 
349.8 
266. 1 


280.2 


ፎን 


AO 
416.8 


299.7 


230.1 


147.7 
309.0 
331.7 


d 


1595 
181.4 


194.1 


208.6 
142.2 
lon 
118.6 
15728 


150.3 
1817 
193.3 
202.8 
141.9 
130.0 
PS 
161.6 
182.8 


Po/Pa 

7 

20927 
1.996 
1.997 
1993 
2.005 
1.974 
1.990 
28977 
2.006 
23. 000 


= 


790 
916 
, 70۵ 
.941 
963 
772 
8 
ናኝ 


S‏ م 
ው C‏ 


e pe pr po pp. PJ] Pp نم زح عم هم عمو عن‎ 
to 


ot 


ZI 


2.961 
2.940 
1,957 
1.962 
1.960 
Dy» 


3.994 


VPRESS 
Vrms 
0.28391 
0.18870 
0.18890 
0.18860 
0.18940 


18982 
2835 

18873 
18905 
18947 
18939 
283540 
28567 
18900 
09465 
18886 
18922 
0. 18899 


كه c9‏ و و و و وه و و و Co gm‏ 


0. 28390 
0. 28310 
0. 28510 
0. 28380 
0. 18870 
0. 28530 
0. 19000 
0.18938 


0. 28380 
0. 28540 
0. 28291 
0. 28291 
0.18710 
0. 18870 
0.18710 
0. 28340 
0. 28421 


PMEAN 
psia 
128, 26 
149.97 
150, 00 
149.96 
149.97 
1518597 
150. 50 
149.90 
149.94 
150. 08 


149, 31 
150. 30 
126, ረፅ 
ን 
154, 80 
133, 00 
150. 15 
150. 05 
149.96 
120.14 
120, 88 
150.17 
120, 11 


148.01 
147. 20 
146. 76 
146. 31 
146.91 
148.26 
148. 71 
147.90 


191.49 
151. 68 
۳ 
152. 60 
152.43 
SA 
153, 03 
120, 24 
124.40 
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PMEAN 

aVrms 
100. 17 
99898 
100. 00 
29.97 
99.78 
191.06 
100. 48 
ማገ) 
99.96 
100. 5 


99.54 
100. 0 
104. 17 
102.5 
105. 20 
102.00 
180, 10 
100. 04 
99.97 
100. 07 
180, 03 
100.11 
100. 07 


98.67 
98.13 
77.84 
37.94 
37.94 
98.84 
33.14 
38. 60 


100.99 
101. 12 
181, 83 
101.73 
101.65 
181, 71 
102.03 
100. 16 
89.60 


FREO 
Hz 
22.05 

923. 20 
37.10 

948. 70 

998. 70 

3/0.40 

3/0. 40 

939.50 

318, 10 

208, 40 


307.94 
305.90 
355. 70 
344. 90 
324. 49 
207.32. 
253.28 
301.57 
307.10 
344. 30 
334. 62 
253. 44 
RASO 


239. 30 
239. Bb 
244.21 
248.82 
230.88 
234.76 
258.85 
260. 29 


229.24 
238. 28 
242.48 
246.82 
248.93 
243.35 
234.37 
229.07 
221/21 


TC 
deg € 
"dr 
23/07 
-44.7 
-34.0 
-24.7 
-5 


"5709 
-68.6 
-44.7 
-34,0 
-2423 
zig 
-43.1 
ን 
-68.6 
-34.0 
AT 
-13.4 
-44.6 


TH 
deg C 
لك‎ 
2 
22. 


QA 


PI 
CA 
F3 CO N e wN N AOC س ضط‎ 


۲4 F5 F3 ۰ 


نم 
J O C4 O — O- Cn F‏ 4= 


۲4 ۲.۵ rn 
C4 C4 C4 E ዮጋ cn e La 
ጨ ل‎ 


ኩ.ጋ 


et 
24.1 
24.2 
24,3 
22.4 
2592 
21.0 
LEE 


220 


PD PD PD PO PJ 
PO ل‎ = iN = 
XD owe 


r3 n3 
ied سا‎ 
cl CN -O C4 — cn 


d c 


C 


old 

stack 
pure He 

short 


resonator 


new 
stack 

pure He 
short 


resonator 


old 
stack 
He-Xe 
short 


resonator 


new 
stack 
He-Xe 
short 


resonator 


HEATLK 


watt 

0.3692 
0,5413 
0, 4576 
0.3869 
0.3297 
0. 2298 
0. 2482 
0. 4740 
0.5780 
0.5474 


9.5075 
0.3934 
0.4278 
0.3623 
9. 5056 
0.2369 
0.4410 
0. 522: 
0.5746 
۵, 2445 
0. 5057 


7113 
9 à» لا‎ 


۵, 4297 


0. 5226 
0. 5414 
2.4870 
0. 4288 
0.5271 
0.5971 
d. 5544 
0, 2241 


8.5121 
d. 5161 
8. 4624 
3.4105 
9,654 
9,4410 
ዐ, ጋዓኋ። 
0.5177 
0.5320 


COPR 


ci7/carn 


0. 0000 
0. 0200 
0. 0511 
0.0475 
à. 0551 
0. 047? 
0. 0506 
0.0380 
0. 0800 
0. 0300 


0, 0000 
0, 0000 
0. 0481 
0. 0508 
0.0612 
0.0398 
0.0635 
û. dd 
à. 0000 
0. 0000 
0.06721 
0.055 

0.0556 


0. 0000 
à. 0404 
0, 0620 
0.0767 
0.0752 
0, 0000 
à. 0000 
0.0797 


9. 8088 
0. 0587 
9. 0856 
0. 0983 
0.0875 
0,009: 
Û. 0000 
9. 0000 
0. bude 


CUP 


Q=heater 


0. 0000 
0, 0000 
ሀ, 1033 
0.1998 
0. 2860 
D. 3629 
0.3379 
d. 1714 
0. 0000 
0. 0090 


9. 0000 
0. 0000 
0.1619 
0. 2927 
8.3185 
0. 3072 
0. 2088 
0, 0000 
0. 0000 
0, 0000 
0.52 

8, 2837 
0. 1688 


. 0000 
1013 
1871 
2638 
2902 
0006 
. 0000 
SS 


ር ፻ ርጋ ርመ 0 ር3 ርጋ ሮ። 


ሀ, 2000 
0.1727 
0.529 
0.3501 
0. 3407 
0.2139 
0. 0000 
۵. ۵ 
۵, 00 


WORK 

watt 

259719 
4,5736 
4, 7429 
4,9991 
J. 2348 
3, 3306 
J. 97/47 
18. 5129 
4,1944 
9:997) 


4. 2636 
9. 2190 
6677 
1953 
2307 
. 8244 
. 2883 
SO 
.2091 
. 3606 
. 6158 
, 0303 
9701 


pL 
4 Cm CH o e -O — Cn Cn Cn حم‎ 


, 6249 
22077 
22508 
6902 
. 1844 
121/2 
1716 


Ps 


ra rJ 4» ዮየ-3 CN cn عط طط‎ 


ra r3 ra cn en 4> e 
ጾ ፌባ 
ርባ 
O 
co 


9 
4 
A 
C.4 
El 


CHEAT 


watt 

0. 0000 
0. 0000 
0. 5002 
0. 2998 
1.4969 
2.0064 
2.0054 
2 
0.0000 
۵, 0 


9, 0000 
0. 9000 
0.7362 
EI 
67/5 
E 
+ 3060 
9. 2000 
0, 0000 
9, 9000 
1,5006 
2. 5258 


0.2546 


حم — J‏ لم 


0. 0000 
0-4 
0.9772 
1.5010 
0. 7499 
0. 0000 
9, 9000 
1. 4964 


0, 900 
0.3524 
1.5420 
1.8796 
0.3820 
0.5095 
0. 0060 
0. 0000 
0. 0000 


149 


TC/TH 


0.7194 
8. 7307 
9. 7725 
0. 8078 
0.8384 
0. 8848 
0.8768 
0.7576 
0.7123 
0.6703 


0. 6754 
9. 5240 
9. 7709 
0, 5062 
0.8366 
0.8735 
0.7567 
0.5690 
0.6916 
0. 8139 
0.8376 
8.8758 
0.7717 


9, 6851 
9. 7154 
0.7441 
9. 7748 
0.7945 
0. 6832 
0.7161 
0. 8556 


ወ, 5203 
0. 7154 
0. 7448 
0.7739 
0, 7934 


0. 


0.6770 
0.6585 
0. 


ACC 
cgsE-3 


320. 28 
211.14 


GS 2 
O 


Act q 
des نا‎ 


269.80 


ا2 


294.41 


Tey} Q7 


7 لأ ه هل لب 


197.49 


277,10 


197,77 


296.18 
ان 


299. 04 


271.47 


308. 47 


372.86 


P3291 
194,91 
134, 34 
285.27 


33.29 


242.98 


66. 50 
72.68 


ale 


86. 70 


9d 


63.04 


uuo 


67.48 


RHEAT 


ohms 


24, 00 
24, 00 
24.40 
24,48 
94. 55 
74,64 
94.40 
94.00 
74.00 
94, 48 
34. 39 
94, 64 
24. 40 


74. 00 
94.16 
24, 23 
74.41 
94, 45 
74, 32 
24.00 
74, 00 
74. 00 


IHEAT 


mA 
0.00 
0.00 
72.81 
102.82 
129.92 
145. 50 
143.50 
145.00 
0.00 
0. 00 


9, 00 
9. 00 
99. 90 
116.90 
1 3.2 
158. 50 
128, 00 
0, 00 
0. 0 
0. 00 
138. 00 
126, 50 
99.50 


0. 00 
E97 
102.94 
125. 09 
89.11 
0. 00 
0, 00 
Dec 


0 00 
75. 70 
119.39 
141.10 
96.80 
73.30 
9, 90 
۵. 0 
3, 00 


010 
stack 
pure He 
short 
resonator 


new 

stack 
pure He 

short 
resonator 


old 
stack 
He-Xe 
short 
rezonator 


new 
stack 
He-Xe 
short 
resonator 


THERM PEN 


cm 

0. 010133 
0.010085 
0. 009957 
0. 009857 
8.009775 
0. 0095/5 
0, 507621 
0.010005 
0.01014; 
3.01025 


0. 0102627 
0. 0102926 
à. 00273765 
0.009777 
d. 209653 
0.00955 

0. 019038 
0. 010349 
0. 010249 
0.009853 
0, 009314 
0.002717 
9, 087751 


0.012596 
0.01252 

0.012452 
0.012328 
0.012225 
0. 012575 
9.012421 
0.011955 


0.912685 
d. 61237 

0.012272 
0.812134 
0.012040 
0.012177 
3.012274 
0. 012572 
0.013451 


VISC PEN 
cm 
0. 008307 
0. 008258 
0. 038152 
0. 008081 
ዐ, 00501. 
0.007848 
9. 007511 
0. 008203 
0. 005515 
0.028410 


0.008413 
0. 008440 
0. 008031 
0. 602017 
9. 607921 
0.007917 
0. 0038227 
0. 020485 
Ø. 008407 
0. 002077 
d. 008045 
0. 007766 
0. 0081979 


0. 005555 
0. 005498 
0.005450 
0. 085401 
0.005541 
0. 005234 
0.905444 
9. 005202 


à. 006532 
o. 006477 
d. 005557 
0. 005295 
0.006246 
0.005517 
0.005420 
0. 005575 
0. 005778 


150 


DELTA 1 
deg 
85. 70 
79.68 
67.38 
26 70 
47.90 
223800 


ም i= 
E ኤታ 


59.70 
. 
. JÜ 


~O CO 
cn cn 


30. 18 
94.5 
57.38 
57.50 
38.52 
37.50 
76.80 
98.80 
91.20 


c 
34. 7 


46.28 
37. 10 


94.6 


0 
59.5 


75.30 


Bo. 


ناك . 


p= 

i 

>» E 
a 


~ tr E 


-4^ CH -E 
Lei CÙ à 
Cy E 


tot 
watt 
0. 5692 
0.3413 
0.2578 
esos) 


ORME 
ES 


!.፲ ር) E C3 r3 n — سم‎ — A E 
ርት 
E) 
ርር) 
cn 


e AJ 
Cri 
ርት- 
c 
Cr- 


«3226 
0428 
. 4352 
9298 
1391 
fe 
 Jod4 


TCE 
T ga 


p= ኤኤ عسو بم‎ [fX 


—— C Cu 


0.5151 
1.3704 
4.0044 
2, 2982 
1.2555 
0, 7206 
8.5455 
0.5177 
0.62520 


COFRE 
co0/carn 
0.0224 
0.0456 
0.0995 
0.0659 
0. 0571 
0.0526 
0. 0558 
0.0716 
0.0361 
0.0548 


2, ዉ25: 
0. 0237 
0.075 

0.0778 
0,0724 
à. 0527 
0.0734 
0.0335 
93, 9597 
0.0579 
0.0725 
0, 9605 
d. 0518 


0.0331 
0. 0840 


0,0767 . 


8.0986 
0.1141 
0.0377 
00973 
1. 


تع ات 


),ህ/14 
1698 
24189 
120; 
1241 
1295 
. 1140 
O72! 
4.0670 


Gi E 


Cc G& Es E Cz 


ar 
COPE 


w ci? 


0.0574 
l208 
d. 2028 
0. 2772 
d, 3482 
0. 4043 
0, 4 
0. 2355 
0.1291 


6.1331 
0.0623 
ቁ, 2332 
1 

s ند‎ 


ü. 3706 
3. 4478 
9. 2478 
0. 0677 
ሀ, 1563 
0. 25353 
0. 
0. 145 


መ. መ መዋ መና 
Dcus 
Sas whi 


9, 1251 
9, ሪ111 
d. 2813 
doom? 
, 2407 
1305 
2461 


7 alo 
OOS 


Ca Cm Cu E 


ር) 
E. 
[3 
CI 


ር) 
[A4 1 
cr 
p 


th a 
-. a لي‎ . a 


old 

stack 
pure He 

short 
resonator 


ner 
stack 

pure He 
short 

resonator 


old 
stack 
He-£e 
short 
resonator 


በ=8 
zberi 
SVEN 
۳۵-۸ 
short 

resonator 


VHEAT 
V 


VLOCKIN 
aVrms 


198. 
208. 
ረረ 
20d. 
160. 
147, 
153. 
| 


“ያዛ - CO ርኩ ርር) Cd  ጐኤ 


212. 
249, 
275, 
298 " 
326. 
212. 
206, 
425. 
425. 
316. 


h3 CA —J —1 t9 cO 


~J e Cr E) 


183. 
G: 
208. 
226. 
137, 
142.0 
126. 
182. 


ሯ- -4 y un 


— pa 


Po/Pn 


Í 

2. 968 
2.778 
2. 968 
24157 
12997 
no 
235 
2. 988 


2772 
por 
. 001 
. 003 
994 


4 r3 PI ዮጋ ሦ* PJ PJ) € ሁሙ 


led PA PO PO na عو هم زيط هو هو‎ 
~O 
co 
0 


VPRESS 
Vrms 
ፅ, 28511 
ሰ, 28201 
0. 283521 
0.28303 
0. 18950 
0.18878 
0.18921 


0. 28523 


0. 18715 
0. 18902 
0. 18920 
0. 18902 
0.18954 
8.189352 
0. 28343 
0. 23363 
ወ, 28357 
0. 28335 


0. 28315 

. 28300 
0. 28310 
0. 28258 
0. 18742 
«15225 
Jig 
29274 


© لت بهت 


. 0 
. 18702 
. 18948 
. 18915 
18957 
. 18923 
. 18241 


go ሮ3 GO O ሮ3 ሮ3 C» 


02 
0. 13858 
0. 18978 
0. 198856 
0. 18700 
0. 2832 

0. 28546 
0. 2869 
0. 28529 


8. 17804 


PMEAN 


psia 

131, 26 
138, 74 
51 2ت‎ 
13ረ, 76 
19125 
130. 38 
158. 38 
130. 33 


150. 60 
149.97 
149.93 
149,78 
150.57 
138, 13 
149. 88 
120, 80 
149. 90 
149.97 


149.91 
149.73 
149.55 
149.46 
120, 14 
149.99 
150. 03 
149.97 


149.67 
149.91 
149.87 
147, 72 
149,87 
150. 02 
149. 78 


120, 82 
147, 76 
149.73 
150. 84 
IS ند‎ 
1350. 41 
150. 24 
130, 27 
150. 30 
00002 


Pn 


PMEAN 
aVrms 


100. 4 
100. 49 
108. 99 
101.84 
100. 82 
100. 25 
108. 25 
100. 25 


120, 40 
29.99 
du 
99.80 
100. 8 
100. 10 
ከ ገገ. 
100. 00 
Weds 
9008 


79, 94 
22.82 
99.7 
99.54 
100. 07 
50 
108. 02 
pod 


8 
994 
57871 
Duro 
901 
100.01 
77... 


100. 01 
79.04 
9 
100. ጋዕ 
100.35 
100, 27 
18, 16 
100, 15 
100. 20 
39. 68 


FRED 
Hz 
ድቃ. 
299,42 
305.07 
318.85 
313. 60 
306.06 
296. 28 
291.54 


470. 43 
483, 36 
495. 98 
304. 18 
314. 53 
459. 90 
488. 19 


د ا 
99.5 


498. 52 
458.88 


294. 28 
206.60 
201171 
276.91 
ሬን 
20 00 
286. 19 
2835.80 


20797 
21:599 
219777 
222.90 
ያር 
231955 
207.98 


490.81 
5 ../ 
386, 27 
SE 
ሣራ? 
317.66 
307. 60 
498.35 
493.36 
472. 10 


21209 
-60. 4 
-51.8 
-42.9 
2 
-00. 1 
-65.0 
x 


1 | | 
ci ظط‎ O 
J> & E) 

CA 


1 
3 
cn 


1 
عم‎ 
ted 


۱ ۱ 
C-i = 
PO ted 


4۸ 
۲۳۰. 


۱ 
0 
ኑጫ የ 
O C4 C4 — ኢ3 -4 - E) CO cC 


BOL. 


| 
p 
eye 
ON O سر خط‎ C4 -J € يع حم ل‎ 


EU 


نم لم ټم h3‏ 3 


F3 = 
] Pa cn 

a a . 9 

GS ro UN CA El rá 


24, 
24. 
29, 
ea. 
ማን 


do d 8 


77) 


y ۵ 
D^ 


P 
£ te 


4 o QQ O ርስ አ پډ‎ ጭ 


zi 


22, 


139 


AY 


ማጣ 
Pau 


p 


2] 


d e 


“Y 0۵ Co hJ - c co 


Z5: 


ZEN 


ata! 


is 


24. 
24. 


ማኘ 


e 


ን 
ፌ 


527 


ais 


26.0 
24,8 


A N GE)‏ هم نم 


م يس هم 


new 
stack 
He-Ar 
(225) 
short 

resonator 


new 

stack 
pure He 

long 
resonator 


new 
stack - 
He-Ar 
(19%) 
long 
resonator 


new 
stack 
He-Xe 
long 

resonator 


dual 

stack 
pure He 

long 
resonator 


HEATLE 


watt 

0. 5200 
0.4718 
0. 4408 
9. 39:8 
0. 3509 
0.42 

0.5059 
0. 0 


سم و صرح 
0.54247 


0.4329 
d. 2538 
0. 5131 
0.2432 

| ES 
d. 4448 
0. 3349 
0.3349 


9, 3339 


à. 3746 
0.2225 
Q. 2005 
0.4477 
0. 5985 
0.4743 
0.5395 


9. 5605 


0. 3-549 
2, 4350 
0.918 
0.3112 


1 2 
ددع وه 


0.1644 


0 CTT 
8 لال ءل‎ 


ED‏ میا 
3776 


8, 
8.3410 
እ. 


COPR 


Q=heater cl7/carn 


Q. 0000 
9, 9422 
0. 0521 
0.0729 
0. 0685 
0. 0505 
۵۰ 00 
9, 0008 


à. 0080 
9. 0898 
0.1145 
0.1172 
0.1071 
0. 000% 
0. 1002 
0.1172 
0.1176 
۵, 0 


0. 0000 
0.035 

0.3980 
0. 1293 
0. 1241 
0.0787 
0. 0000 
a. 0000 


0. 8800 
0. 1871 
0. 1540 
1 E 
0.1405 
9. 1077 
0, 0300 


9. 0008 
0.0317 
6, 0384 
d. 2078 
ó. 0614 
0.0633 
0. 0477 
0.01927 
3, 0000 
8, 9090 


COP 


0. 0000 
0. 1059 
0.1910 
0, 2472 
0. 2637 
0.1346 
0. 0000 
9. 0000 


9. 0000 
d. 2205 
0. 4778 
0. 5851 
9. 7198 
0.0000 
0.3266 
0,4627 
0.4636 
0, 0080 


۵۰ 2028 
8. 0327 
0. 2817 
9. 4525 
0.4731 
0.2412 
0.0000 
0, 0000 


2, 0000 
0.5286 
0.5770 
0.7962 
0. 9840 
1.1182 
0. 0308 


û. 9000 
d. 1512 
3. 2798 
2.3741 
20172 
0, 4021 
3, 2414 
0. 0787 
3. 0060 
9. 2800 


WORK 

Watt 

de 2 182 
3.7066 
6.0427 
6.2969 
2,9544 
2.7462 
2, 4825 


ES 


. 1 


4, 5328 


3. 0525 
J. 1409 
3. 9424 
5,8926 
10. 4281 


11. 5415 
11.5549 


8.6416 


, 0588 
ፈያ 
. 7004 
2051 
9346 
6677 
3261 
1180 


۲.۵ ۲.۵ FJ O- CN CA cn 


محل 


ግማ1 
ፀ/ፌ/ልፈ 


746 
2032 


545 


«0363 
6300 
7336 


C^ F3 r2 r9 m m‏ سر 


E 
ር 63565 


3. 9023 
4,1467 
4,3714 
4,/417 
18. 0371 
۵ 
8.7953 
g. 4337 
5. 0281 


Ld 


GHEAT 


watt 
0. 0000 
0.5041 
1.075 
1.3065 
0.7836 
0. 4247 
0. 0000 
d. 0000 


0. 0000 
1,3394 
2. 5885 
ome! 
4.2776 
0. 0080 
2, ዓረ39 
35242 
3, 3368 
0. 0000 


0,000 
9, 4/20 
1. 6082 
2. 6834 
1.3941 
0.6275 
0. 0000 
۵, 0000 


0, 0200 
ል, 633 

1,642 
0263 
./910 
4459 
. 0000 


cu Crd ra r3 


0000 . 
2922 
2577 
ا 
7259 
0363 


EL) 


0.6738 
. 0060 
0.0000 


ret 4 p ft ኡሙ ES CD 


p 
r 
ር 
ር 
F 


5 2 


TC/TH 


8.575 
9, /128 
0.7444 
0.7723 
0.73592 
Q. 7544 
0. 7048 


0. 6783 


9. 7190 
0. 7687 
0. 8051 
0. 8331 
0.8709 
0,7159 
0.7632 
m7 
0.7977 
0. 7056 


0.7029 


9.7146- 


0.7419 
0.7679 
0.7729 
0, 7354 
9.7189 
d. 6999 


0. 7122 
8. 7542 
0. 7950 
0.8333 
0.8751 
0.9121 
0,7127 


9. 7829 
2,852 


ori 


d. 826 
d. 3743 
d. 7095 
0.8552 
0.8292 
0.7985 
0. /941 
0. 7:05 


ACC 

cgsE-3 
2l 
104. 08 


112226 
lioe 
84.15 


76. 59 


66.36 


OF 7 
2 


Dy «^J 


e 


0 


229.24 


nes E 
ፌዴኒ/4 8 ኒህ 


280. 50 
166.56 
314.77 
394.32 


354.90 


248. 04 


97.14 


21597 
101, 40 
112,45 


19.86 


569.40 
60. 27 


86. 33 


32.68 
5 
33 
31.60 
9 7 


54,72 


326 


162, 75 
175. 54 
125073 
28. 2o 


255 Qt 
ረ ን 


Sales 


292,77 
268.71 


sec 4 
sis e 


278.94 


RHEAT 


ohms 

94. 00 
94.16 
74.27 
94,41 
24.45 
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